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photographed  under  a SEM  and  digitally  measured  according  to  the  Richardson  method 
to  obtain  fractal  dimensions.  The  results  illustrate  a decrease  in  fractal  dimension  with  an 
increase  in  fracture  toughness  for  ductile  materials.  The  data  are  compared  against  two 
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CHAPTER  1 
INTRODUCTION 


Since  the  inception  of  the  science  of  fractography,  fracture  surfaces  have  been  an 
invaluable  source  of  information  regarding  the  fracture  event.  Investigations  into  the 
various  modes  of  failure  have  produced  significant  advances  in  the  qualitative  and 
quantitative  understanding  of  how  and  why  fractures  develop.1"  ’ The  application  of  this 
knowledge  has  increased  the  safety  and  efficiency  with  which  new  structures  and 
mechanical  systems  can  be  created  and  used.  New  investigative  methods  and  techniques 
promise  to  add  more  information  to  this  important  field  of  materials  science. 

In  recent  years  a great  deal  of  interest  has  developed  in  the  area  of  quantitative 
fractography,  that  is,  measuring  the  size  of  fracture  features  or  material  properties  directly 
from  a fracture  surface.  One  of  the  concepts  that  has  emerged  as  a potentially  useful  tool 
in  these  efforts  is  that  of  fractal  geometry. 

Fractal  geometric  analysis,  when  applied  to  a fracture  surface,  provides  a measure 
of  its  irregularity.  To  date,  a large  number  of  studies  using  these  techniques  have  been 
performed  on  a variety  of  fracture  surfaces.  A subject  of  particular  interest  has  been  the 
potential  relationship  that  exists  between  plane-strain  fracture  toughness,  usually 
expressed  as  K.ic,  and  the  degree  of  fracture  surface  tortuosity.  Results  have  been  mixed. 
Some  studies  have  shown  a positive  correlation,  some  have  shown  a negative  correlation 
and  others  still  report  no  correlation  at  all.4 
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A review  of  the  literature  shows  that  the  efforts  made  to  date  in  correlating  K|C  to 
the  fractal  dimension  have  employed  a wide  variety  of  materials,  fracture  modes,  and 
measurement  methods.  In  most  of  these  studies  the  material  is  unique,  the  fracture 
process  is  not  well  understood  and  the  preferred  measurement  method  is  debatable.  The 
variety  and  complexity  of  the  information  obtained  have  made  the  assessment  of  material 
properties-surface  feature  relationships,  and  the  drawing  of  general  conclusions,  a 
difficult  task. 

Despite  inconsistencies  in  the  literature,  some  general  consensus  is  now  emerging 
among  researchers.4  The  separate  elements  employed  in  the  construction  of  a fractal 
dimension-fracture  toughness  relationship  are  becoming  more  formalized  and  accepted  as 
new  studies  demonstrate  their  validity.  The  present  study  is  intended  to  become  part  of 
one  of  the  building  blocks  which  demonstrates  how  fractal  analysis  might  be  applied  and 
what  material  property  information  might  be  gained  from  fracture  surfaces.  In  particular, 
this  study  focuses  on  the  fractal  dimension-fracture  toughness  relationship  in  heat  treated 
AISI  4340  steels.  The  significance  of  any  relationships  that  emerge  will  also  be 
considered. 

General  Comments  on  Fracture  Toughness 

It  is  the  principal  goal  of  every  structural  or  mechanical  engineer  to  prevent  the 
sudden  failure  of  their  particular  component  or  assembly  by  design.  The  consequences  of 
ignoring  the  conditions  that  lead  to  fracture  are  disastrous.  It  is,  therefore,  of  extreme 
interest  to  understand  the  factors  that  may  control  and  predict  this  potential  outcome. 

Unstable  fracture  in  metals  typically  begins  from  some  preexisting  flaw  or  crack. 
Cracks  may  develop  in  a component  for  a number  of  reasons,  fatigue,  stress-corrosion 
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and  hydrogen  embrittlement  being  but  a few.  Regardless  of  how  the  flaw  initiates  and 
grows,  unstable  crack  propagation  is,  in  general,  the  final  event  that  causes  structural 
failure.  A knowledge  of  when  a component  is  approaching  this  end  condition  is  vital  in 
original  design  and  failure  prevention.  Furthermore,  a knowledge  of  what  material 
features  assist  in  delaying  or  preventing  this  occurrence  is,  for  obvious  reasons,  crucial. 

The  parameter  most  commonly  used  to  describe  a material’s  resistance  to  unstable 
fracture  is  fracture  toughness.  It  is  generally  represented  by  values  of  the  plane-strain 
fracture  toughness  factor  KiC.  In  practice,  a structure  containing  a flaw  can  assume  any 
value  of  stress  intensity  at  a crack  tip  up  to  this  critical  quantity.  When  it  is  reached,  the 
component  fails  catastrophically.  Stress  intensity  depends  on  applied  stresses,  flaw  size 
and  flaw  geometry.  The  critical  stress  intensity  is,  in  general,  a material-related  property. 

The  value  of  fracture  toughness  (critical  stress  intensity  factor,  Kic)  for  a 
particular  material  is  commonly  found  through  standardized  testing  methods.  It  is  now 
well  known  that  a number  of  geometric  and  microstructural  features  such  as  specimen 
size,  loading  geometry,  crack  geometry  and  orientation,  grain  size,  second  phase  particles 
etc.,  will  affect  the  measurement  and  apparent  value  of  this  parameter.  It  is,  therefore,  of 
practical  and  scientific  interest  to  determine  how  material  characteristics  influence 
fracture  toughness,  the  various  methods  through  which  this  property  may  be  properly 
measured  and  how  the  fracture  mechanism  itself  is  related  to  the  critical  stress  intensity 
factor.5 

Significance  of  Fracture  Surfaces 

Just  as  research  on  the  relationship  between  fracture  toughness  and  microstructure 
has  added  new  knowledge  about  the  mechanism  of  fracture,  investigations  into  the 


4 


relationship  between  fracture  surface  features  and  fracture  toughness  offer  to  do  the 
same.  There  is  reason  to  think  that  these  two  factors  are  related.6 

An  examination  and  comparison  of  metallic  fracture  surfaces  created  by  overload 
reveal  that  there  is  an  apparent  connection  between  fracture  surface  roughness  and  the 
associated  fracture  toughness.  That  is,  materials  with  relatively  smooth  fracture  surfaces 
appear  to  have  a lower  resistance  to  fracture.  This  notion  seems  to  be  intuitively  correct. 
It  appears  logical  for  relatively  smooth  surfaces  to  consume  little  energy  in  their 
formation  and,  thus,  reflect  this  fact  in  low  values  of  fracture  toughness.  The 
mathematical  tools  are  now  available  to  determine,  in  a quantitative  way,  if  this  notion  is 
indeed  correct  and  what  that  relationship  might  be.  The  connection  between 
microstructure,  fracture  surface  texture  and  fracture  toughness  is  expected  to  lead  to  a 
better  understanding  of  the  fracture  process,  the  development  of  tougher  alloys,  and  the 
identification  of  pre-fracture  material  properties  after  a component  has  failed. 

In  order  to  study  the  relationship  between  fracture  toughness  and  the  resulting 
fracture  surface  topography,  fractal  geometry  will  be  used.  Fractal  geometry  is  a non- 
Euclidean  geometry  that  can  be  employed  to  quantitatively  describe  the  tortuosity  of 
fracture  surfaces,  which  are  characterized  by  their  fractal  dimension.  A 4340  steel  will  be 
heat  treated  to  vary  the  fracture  toughness  and  correspondingly  change  the  fracture 
surface  topography.  These  heat  treatments  will  cause  the  behavior  to  vary  from  fracture 
in  a brittle  manner  to  fracture  in  a ductile  mode.  The  results  will  be  compared  to  existing 
and  new  theories  to  explain  the  observed  relationship  between  toughness  and  the  fractal 


dimension  of  the  fracture  surface. 


CHAPTER  2 
LITERATURE  REVIEW 

Fractal  Geometry 

Concepts 

The  concept  of  fractal  geometry  is  relatively  young.  Fractal  geometry  expands 
the  concept  of  dimension  and  recognizes  that  there  are  an  infinite  number  of  dimensions 
in  between  the  topological  dimensions  of  1,  2 and  3.  In  the  particular  case  of  irregular 
surfaces,  for  instance,  it  is  possible  to  obtain  fractal  dimension  values  between  2.0  and 
3.0.  This  number,  to  a large  extent,  represents  the  degree  of  surface  irregularity.  A 
planar  surface  would  be  expected  to  produce  values  near  2.0  while  more  tortuous  surfaces 
would  be  expected  to  be  closer  to  3.0. 7 

Although  rougher  surfaces  are  thought  to  exhibit  a greater  fractal  dimension,  this 
statement  is  not  strictly  true.  The  industry  standards  that  cover  roughness  usually  define 
it  as  a number  which  represents  an  average  deviation  from  a mean  line.  For  example;  the 
International  Standards  Organization  (ISO  4287)  defines  the  most  commonly  used 
roughness  parameter  Ra  (roughness  average)  as8 

L 

(2-1)  Ra=l/L  j|  y(jc)  | <&■ 

o 

This  expression  is  illustrated  graphically  in  Figure  2-1. 
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Figure  2-1 . Representation  of  roughness  average  (Ra)  per  ISO  4287. 
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Figure  2-2.  Lines  of  similar  Ra  but  different  fractal  dimension. 


Consider  the  two  lines  in  Figure  2-2.  If  the  traditional  engineering  definition  of 

roughness  is  applied  to  these  lines  then  Ra  would  be  found  to  be  very  similar.  Both 

profiles  deviate  about  the  same  amount  from  a mean  line.  If  their  fractal  dimensions 

were  measured,  however,  these  values  would  be  expected  to  differ  significantly.  Line  2 

is  considerably  more  tortuous  (irregular)  than  line  1. 

An  accurate  definition  of  fractal  dimension,  then,  is  not  based  strictly  on 

roughness.  Fractal  objects  possess  at  least  these  two  important  qualities: 

1 ) Scale  invariant  self-similarity:  Fractal  objects  display  identical  levels  of 
irregularity  no  matter  what  the  magnification  of  the  object.  In  fact,  given  two 
images  of  the  same  object  at  vastly  different  magnifications,  the  viewer  would  not 
be  able  to  determine  the  actual  scale  of  the  object  without  a reference  scale. 
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2)  The  geometric  features  of  length,  area  or  volume  change  at  a predictable  linear 
rate  on  a log-log  graph  when  the  measurement  device  changes.  For  example,  if  an 
irregular  line  is  measured  with  a step  length  of  1 cm  (using  dividers,  for  instance) 
the  total  length  will  be  longer  than  if  a step  length  of  2 cm  is  used. 

Measurement  Through  Richardson  Plots 

The  science  of  measuring  fractal  curves  made  an  important  advance  when 
Richardson  attempted  to  correlate  the  length  of  national  boundaries  with  their  respective 
military  and  economic  conflicts.  The  boundary  measurements  were  made  from  various 
maps;  however,  he  quickly  realized  that  these  perimeter  values  varied  widely  depending 
on  the  scale  of  the  map.  He  further  learned  that  boundary  dimensions  also  varied  if  the 
scale  was  kept  constant  and  the  measuring  unit  was  altered. g 

Figure  2-3  shows  Richardson’s  measurement  process.  Each  of  the  four  identical 
images  produces  a different  boundary  length  depending  on  the  ruler  length. 


Figure  2-3.  Variation  of  boundary  length  with  ruler  length.  Smaller  rulers  produce 
larger  perimeters  [Ref.  4,  pg  27,  with  kind  permission  of  Springer  Science  & 
Business  Media]. 
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The  variation  of  perimeter  with  measurement  scale  can  be  demonstrated  with  a 
log-log  graph.  The  typical  result  is  a linear  relationship  as  shown  in  Figure  2-4.  Plots  of 
this  type  are  often  used  today  in  the  fractal  analysis  of  irregular  lines  or  perimeters  and 
are  referred  to  as  Richardson  plots. 


Figure  2-4.  Richardson  plot  produces  a straight  line  on  a log-log 
graph.  The  slope  gives  the  fractal  dimension.  Note  that  the 
deviations  from  linearity  at  the  ends  are  due  to  measuring  unit  effects. 

[Ref  4,  pg  29,  with  kind  permission  of  Springer  Science  & Business 
Media], 

The  mathematical  representation  of  a Richardson  graph  (for  a boundary  line),  and 
the  one  which  is  characteristic  of  fractal  relationships,  is  given  by 

(2-2)  L = Ks(I‘D) 

The  linearized  version  is 

(2-3)  Log(L)  = Log  (K)  + (l-D)Log(s) 

where  in  both  (2-2)  and  (2-3) 


L = Measured  length. 

S = The  measurement  scale. 


K = A constant. 

D = The  fractal  dimension. 
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Other  Measurement  Techniques 
Dimensional  analysis 

Mandelbrot  and  co-workers 10  were  the  first  to  examine  fracture  surfaces  using 
techniques  derived  from  Richardson’s  observations.  For  the  case  of  fracture  surfaces,  the 
Richardson  relationship  was  modified  to  read  as  follows 

(2-4)  A = Kp2/D 

The  linearized  version  is 


(2-5)  Log(A)  = Log(K)  + 2/Ds  Log(p) 

where  in  both  (2-4)  and  (2-5) 

A = Measured  area.  K = A constant, 

p = The  perimeter  length.  Ds  = The  surface  fractal  dimension. 

2/Ds  = The  slope  obtained  from  the  Log(A)  - Log(p)  plot. 

In  practice,  this  technique  is  applied  by  mounting  and  polishing  a fracture  surface 
in  such  a way  as  to  allow  only  the  cross  sections  of  "high  spots"  or  "islands"  of  the 
fracture  surface  to  show  on  a sectioning  plane.  Figure  2-5  illustrates  an  example  of  a 
group  of  such  islands. 


Figure  2-5.  Typical  cross 
section  through  a fracture 
surface  revealing  several  islands 
and  low  spots. 
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Once  a set  of  islands  is  obtained,  the  surface  may  be  digitized  and  a computer 
program  may  be  applied  to  obtain  data  measuring  areas  and  perimeters  for  all  the  islands 
available.  The  data  is  then  used  to  obtain  graphs  of  log  area  v log  perimeter. 

The  Minkowsky  dimension 

The  Minkowsky1 1 dimension  is  determined  for  a line  or  boundary  by  sweeping 
out  the  feature  from  beginning  to  end  with  circles  of  various  radii.  A graph  of  log  area  of 
the  circle  v log  radii  of  the  circle  produces  a slope  which  gives  the  Minkowsky  fractal 
dimension. 

The  Kolmogorov  dimension 

The  Kolmogorov  dimension  is  determined  for  a line  or  boundary  by  covering 
the  feature  with  grids  of  different  sizes.  A graph  of  log  grid  size  v log  number  of  grids 
through  which  the  feature  passes  produces  a slope  which  gives  the  Kolmogorov  fractal 
dimension. 

Fourier  analysis 

A log-log  graph  is  constructed  of  the  magnitude  v.  frequency  for  the  Fourier 
transform.13  This  graph  produces  a straight  line  whose  slope  is  related  to  the  fractal 
dimension. 

Selection  of  Measurement  Technique  and  Inherent  Differences. 

The  fractal  dimension  measurement  techniques  listed  above  are  not  the  only 
methods  available;  however,  they  are  the  most  widely  used.  Despite  the  fact  that  they  all 
claim  to  measure  a "fractal  dimension"  by  definition,  the  only  method  that  produces  a 
true  fractal  dimension  (Hausdorf  dimension)  is  the  Richardson  plot.4  The  other  methods 
have  been  developed  and  employed  in  cases  where  the  feature(s)  are  more  easily 


measured  using  one  of  the  alternative  techniques.  In  addition,  in  many  cases  the  exact 
value  of  the  fractal  dimension  is  not  necessary.  The  important  aspect  being  studied  is  the 
difference  in  fractal  dimension  between  similar  features  or  their  rate  of  increase.  As  long 
as  the  measurement  technique  is  consistent  among  the  features  examined,  and  its 
strengths  and  limitations  are  understood,  any  of  these  techniques  may  be  used.  An 
excellent  review  of  the  available  methods  and  their  inherent  strengths  and  weaknesses  is 
provided  by  Russ.4 
Self-Similar  v.  Self-Affine  Surfaces 

One  of  the  problems  encountered  in  the  literature  is  the  erroneous  application  of 
fractal  dimension  measurement  techniques.  It  is  common  for  many  authors  to  use  profile 
sections  from  fracture  surfaces,  apply  the  Richardson  technique  to  the  resulting  line  and 
report  a fractal  dimension  - property  relationship.  The  error  here  is  that  there  is  no 
evidence  to  suggest  that  the  fractal  dimension  of  a profile  is  equivalent  to  the  fractal 
dimension  of  the  surface. 

The  argument  initiated  above  is  particularly  troublesome  in  the  case  of  fracture 
surfaces.  These  surfaces  are  thought  to  be  self-affine  rather  than  self-similar.4-14  That  is, 
a self  affine  surface  exhibits  different  scaling  behavior  with  orientation.  The  patterns 
formed  on  the  X-Y  (horizontal)  plane  are  not  thought  to  be  the  same  as  those  formed  on 
the  X-Z  (vertical)  plane.  Even  though  some  relationship  may  exist,  they  are  not  thought 
to  be  identical.  In  the  case  of  fracture  surfaces,  the  difference  in  scaling  behavior  is 
attributed  to  differences  in  applied  stresses,  constraints,  residual  stresses  and,  perhaps,  the 


characteristics  of  the  material  itself. 
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If  we  are  to  apply  the  Richardson  technique  we  must  use  a plane  that  is  essentially 
parallel  to  the  direction  of  fracture.  Such  planes  are  known  as  zeroset  planes  and  are 
considered  to  be  suitable  for  a fractal  dimension  measurement  which  is  representative  of 
the  fracture  surface.  Deviations  of  5°  or  more  from  this  zeroset  plane  may  lead  to  large 
measurement  errors. 1 5 

Application  of  Fractals  to  Material  Surfaces 

Fractal  Mechanics 

The  previous  section  provided  an  introduction  to  the  concept  and  mathematics 
supporting  fractal  dimension  measurement.  A number  of  excellent  references  are 
available  in  the  literature  which  describe  the  methods,  underlying  theory  and  applications 
in  greater  detail.  Considering  that  a thorough  review  of  the  underlying  mathematics  is 
beyond  the  scope  of  the  present  study,  the  reader  is  referred  to  those  papers  for  additional 
information. 16-25  This  document  will  elaborate  in  greater  detail  on  those  papers  and 
references  which  are  most  relevant  to  the  questions  evaluated  by  this  work. 

Many  of  the  references  cited  in  this  study  employ  one  the  two  forms  of  the  slit 
island  analysis  method  (SIA  by  Richardson  plots  or  dimensional  analysis,  a.k.a.  the 
perimeter-area  method).  A thorough  review  of  the  mathematical  foundations  of  these 
methods  and  their  application  to  fracture  surfaces  is  provided  by  Meisel.26  Meisel  also 
applies  these  methods  to  "artificial"  fractal  objects  and  determines  their  validity,  strengths 
and  weaknesses.  The  central  finding  in  Meisel's  work  is  that  perimeter-yardstick 
calculations  are  essentially  equivalent  to  perimeter-area  calculations  provided  that  the 
measuring  units  are  chosen  appropriately.  Under  carefully  chosen  conditions  both  will 
provide  acceptable  measures  of  the  fractal  dimension.  Meisel  offers  guidelines  on  the 
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selection  of  length  scales  for  obtaining  satisfactory  measurements.  In  particular  he 
suggests  that  a fairly  large  range  of  ruler  lengths  be  employed  for  measurements  such  that 
scale  dependent  behavior  may  be  detected  in  the  fractal  dimension  graphs. 

Application  of  Fractals  to  Ceramics 

A significant  portion  of  the  work  performed  on  fracture  surfaces  has  been 
accomplished  on  ceramics.27'28  For  this  class  of  materials,  much  of  the  research  has  been 
performed  by  Passoja,  Mecholsky  and  co  workers.29'27  It  is  of  particular  interest  that 
these  studies  did  not  simply  limit  themselves  to  identifying  that  a surface  is  fractal  but 
continue  on  to  suggest  a structure-property  relationships  to  explain  this  behavior.  The 
fundamental  relationship  suggested  by  these  researchers  is  that  fracture  toughness  may  be 
related  to  the  fracture  surface  topography  as  follows 

(2-6)  KIC  = A(D*)I/2 

where 

KIC  = Plane  strain  fracture  toughness. 

A = a parameter  characteristic  of  the  material  class. 

D*  = the  fractional  portion  of  the  fractal  dimension. 

Mecholsky  et  al.30  proposed  that  A is  a function  of  E,  Young's  Modulus,  and  ao,  is 

a characteristic  length  such  that  A = E(ao)1  2.  Mecholsky  et  al.31  note  that  other  functions 

may  fit  the  experimental  data,  however,  equation  2-6  was  constructed  on  the  basis  of 

dimensional  arguments.  The  resulting  units  are  MPa  root-meter  as  required  for  the  plane 

strain  fracture  toughness,  Kic-  The  usual  form  of  their  equation  is  then 

(2-7)  K1C  = E(a<)D*)1/2 

Mecholsky  and  co-workers  have  interpreted  ao  to  be  a length  characteristic  of  the 
material  (a  structure  parameter).  The  actual  meaning  of  this  variable  remains  to  be 


14 


investigated  further.  The  structure  parameter  does  not  seem  to  be  related  to  grain  size, 
distance  between  inclusions,  distance  between  internal  flaws  or  other  typical 
microstructural  features.  The  suggestion  is  made  that  a<>  could  represent  more  esoteric 
"distance  related"  features  such  as  free  volume,  glass-crystal  clusters,  stretched  bonds  or 
glass  phase  stretched  bonds34  in  inorganic  glasses  and  glass  ceramics  (see  Figure  2-6). 


Material  Class 

D* 

Glasses 

0.07-0.1 

Glass  Ceramics 

0.06-0.3 

polycryst.  ceramics 

0.06-0.35 

single  crystals 

0.07-0.12 

ao 

Interpretation* 

10-20 

free  volume 

20-80 

glass-crystal  cluster 

3-20 

stretched  bonds  or 

glass  phase 

1-10 

stretched  bonds 

Figure  2-6.  Relationship  between  D*  and  ao  and  interpretation  of  the  structure  parameter 
for  various  classes  of  materials.  Note:  units  of  a«  in  angstroms  [Ref.  34], 

Despite  the  uncertainty  surrounding  the  interpretation  of  the  structure  parameter, 
the  relationship  (equation  2-7)  appears  to  work  well  for  a large  variety  of  ceramic 
materials.  It  is  particularly  applicable  to  establishing  a relationship  within  classes  of 
ceramic  materials  (see  Figure  2-7). 

In  addition  to  the  fracture  toughness-microstructure  relationships  discussed  above, 
Mecholsky  and  Freiman3"  have  shown  that  there  is  a relationship  between  the  decimal 
portion  of  the  fractal  dimension  (D*)  and  the  flaw-to-mirror  size  ratio  at  the  origin  of 
most  ceramic  fracture  surfaces.  The  relationship  has  been  given  as: 

(2-8)  C/r,  = D* 

A graphical  representation  of  their  data  gathering  procedure  and  findings  is  shown 


in  Figure  2-8. 
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Figure  2-7.  Fracture  toughness-fractal  dimension  relationships  in 
graphical  form  for  various  classes  of  ceramics.  [Refs.  14  and  34], 


Figure  2-8.  Graphical  representation  of  investigative  technique  and 
results  establishing  the  relationship  between  C/ri  and  D*  [Ref  32], 


16 


This  later  finding  is  particularly  interesting  considering  that  several  relationships 
have  been  established  between  mirror,  mist  and  hackle  regions  found  at  the  origin  of 
ceramic  fracture  surfaces/"  Fractal  dimension  studies  imply  that  the  features  that 
develop  at  the  origin  are  somehow  propagated  throughout  the  fracture  surface.  It  is 
assumed  that  all  of  these  features  are  a reflection  of  the  material's  characteristics  and  can 
be  described  by  fractal  analysis  methods. 

Application  of  Fractals  to  Polymers 

The  application  of  fractal  fracture  surface  analysis  to  polymers  and  polymer 
composites  is  more  limited  than  for  metals  and  ceramics.  In  this  particular  class  of 
materials  fractal  concepts  have  found  greater  use  in  the  study  of  polymer  structures  at 
molecular36  and  physical  levels.  7 In  this  respect,  fractal  theory  has  been  applied  to 
polymer  melts, 3X  polymerization  reactions,34  40  rheology41  and  polymer  matrix/second 
phase  composite42-  43  interactions. 

The  focus  on  polymer  structures  not  withstanding,  a number  of  significant  studies 
have  been  performed  which  add  to  our  knowledge  of  the  fracture  process.  Some  have 
emphasized  bulk  mechanical  behavior44’46  while  others  have  delved  deep  into  atomic 
structure-fracture  surface-property  relationships  as  noted  in  the  paragraphs  that  follow. 

Joseph  et  al.47  have  looked  at  the  fracture  surface-fracture  toughness  relationship 
using  atomic  force  microscopy.  Four  relatively  brittle  polymers  were  used  in  this  study. 
Measurements  in  the  mirror,  mist  and  hackle  regions  (identical  to  those  found  in  the 
fracture  initiation  regions  of  ceramics)  resulted  in  fractal  dimension  values  which 
correlated  directly  with  fracture  toughness. 
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Kozlov  et  al.  have  studied  the  applicability  of  fractal  fracture  mechanics  to 
polymers  and  polymer  composites  and  have  concluded  that  the  concepts  should  be 
equally  appropriate  in  this  class  of  materials  as  it  appropriate  in  metals  and  ceramics.  In 
fact,  they  assert  that  the  critical  crack  opening  displacement  method  of  measuring 
fracture  toughness  may  provide  a suitable  scale  of  fracture  in  polymeric  materials. 

Lyu  and  co-workers4  ’ have  applied  fractal  analysis  methods  to  two  engineering 
thermoplastics;  the  partially  crystalline  Polyetherketone-C  (PEK-C)  and  the  amorphous 
polyethersulfone  (PES).  They  report  that  these  polymers  behaved  in  a ductile  manner. 
Presumably,  this  indicates  that  the  materials  fractured  by  crazing  as  is  typical  for 
thermoplastics.  Their  findings  show  a decrease  in  D*  with  increasing  plane-strain 
fracture  toughness. 

Application  of  Fractals  to  Metals 

Fractal  analysis  methods  have  been  applied  to  metallic  fracture  surfaces  from  the 
inception  of  these  concepts.  Mandelbrot,  Passoja  and  co-workers  initially  applied  the 
slit-island  method  and  Fourier  analysis  to  the  fracture  surfaces  of  300  maraging  steel 
separated  by  impact.111  They  observed  that  in  these  ductile  fractures,  the  fractal 
dimension  decreased  with  increasing  toughness  in  a predictable  manner.  This  seminal 
work  launched  fractal  surface  analysis  into  the  mainstream  of  materials  science.  Since 
then,  the  concepts  initiated  by  Mandelbrot  have  spread  to  varying  degrees  to  all  other 
materials  and  their  associated  wide  variety  of  fracture  mechanisms. 

Researchers  in  almost  all  areas  related  to  fracture  have  attempted  to  explain  their 
particular  phenomenon  of  interest  in  terms  of  fractal  concepts.  Fractal  analysis  has  been 
employed  to  explain  the  limit  of  unstable  crack  velocities,50  the  mechanisms  of  fatigue,51 
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55  creep,56  superplasticity,57’58  brittle  fracture,59  and  stress  corrosion  cracking.60  Success 
in  each  of  these  areas  has  been  limited.  None  of  the  studies  have  been  duplicated 
sufficiently  in  other  materials  or  conditions  for  adequate  verification  of  each  theory. 

Pande  and  Richards61  examined  the  fractal  characteristics  of  fractured  surfaces  in 
titanium  using  vertical  SEM  brightness  profiles  and  the  dimensional  analysis  method 
with  slit-islands.  Their  initial  conclusions  indicated  that  the  surface  of  this  material  could 
be  described  using  fractal  geometry.  A second  paper  released  shortly  thereafter  by  Pande 
et  al.6'  cautioned  that  previous  findings  had  been  too  optimistic  and  that  the  fracture 
surfaces  of  titanium  may  not  be  fractal.  In  this  later  paper  Pande  and  co-workers 
attempted  to  relate  the  fractal  dimension  to  the  dynamic  tear  energy  (DTE)  and  found 
only  a tenuous  decreasing  value  of  D*  with  an  increase  in  this  measure  of  toughness. 
Gong  and  Lai63  revisited  this  fracture  toughness-fractal  dimension  relationship  in 
titanium  alloys  by  using  J-R  resistance  curves  and  the  Richardson  method  on  vertical 
profiles.  They  reported  an  increase  in  fractal  dimension  with  toughness. 

A series  of  important  fracture  surface  topography-fracture  toughness 
investigations  have  been  performed  by  Bouchaud  et  al.64, 65  on  ductile  aluminum  alloy 
7475.  In  these  studies,  and  a review  paper  released  later,66  the  authors  argue  that  there 
appears  to  be  a universal  "roughness  exponent"  (£,  related  to  D*  by  D*  = 3 — £ ) which 
measures  in  the  order  of  0.80  ± 0.05.  No  change  was  observed  with  fracture  toughness. 
Other  investigators  have  found  similar  results  in  a variety  of  relatively  brittle  materials 
(plaster,  Bakelite,  porcelain,  graphite,  steel,  and  Al-Si)  which  revealed  a roughness 
exponent  in  the  order  of  0.87  ± 0.07. 67 
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Application  to  Steels 

Relationship  between  impact  toughness  and  fractal  dimension 

Many  researchers  have  long  suspected  a toughness-fracture  topography 
relationship  in  metals.5'6'68  As  noted  earlier,  however,  the  first  attempt  to  identify  the 
relationship  between  toughness  and  "surface  roughness"  through  fractal  geometry  was 
made  by  Mandelbrot  et  al.10  The  study  used  impact  toughness  specimens  of  type  300- 
grade  maraging  steel  and  the  perimeter-area  version  of  the  slit-island  method.  This  study 
was  duplicated  by  Ray  and  Mandal64  (using  a high  strength  low  alloy  steel)  and  Hilders 
and  Pilo70  (using  a ferrite-pearlite  steel)  by  employing  charpy  V-notch  impact  specimens 
and  the  perimeter-area  slit-island  method.  The  findings  of  all  three  studies  are  shown  in 
Figure  2-9. 


Impact  energy  (J) 

Figure  2-9.  Fractal  dimensional  increment  v Charpy  impact  energy  for 
300-grade  maraging  steel  (A),  HSLA  steel  (•)  and  ferrite-pearlite  (o)  steel 
[Refs  10,  69  and  70], 

Hilders  and  Pilo  state  in  their  work  that  their  results  are  consistent  with  the  work 
of  Ray  and  Mandal  yet  opposite  of  that  reported  by  Mandelbrot  et  al.  This  discrepancy  is 
explained  on  the  basis  of  different  micromechanisms  of  fracture.  The  specimens  with 
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increasing  fractal  dimension  as  a function  of  increasing  impact  energy  were  observed  to 
fracture  primarily  through  cleavage.  The  specimens  with  decreasing  fractal  dimension  as 
a function  of  increasing  impact  energy  were  observed  to  fracture  primarily  through 
microvoid  coalescence. 

It  is  interesting  to  note  that  the  results  noted  in  Figure  2-9  are  also  consistent  with 
the  behavior  of  brittle  materials  reported  by  Mecholsky14,34  (see  Figure  2-7;  increasing 
fractal  dimension  with  increasing  plane-strain  fracture  toughness),  consistent  with  studies 
on  brittle  polymer  by  Joseph  et  al.47  (increasing  D*  with  increasing  Kic)  and  ductile 
polymers  by  Lyu  et  al.49  (decreasing  D*  with  increasing  KIC)  and  generally  consistent 
with  the  observations  of  Pande  et  al.62  in  ductile  titanium  alloys  (they  reported  a roughly 
decreasing  fractal  dimension  with  increasing  dynamic  tear  energy-a  measure  of 
toughness). 

Similar  observations  have  been  reported  in  research  where  the  microstructure  was 
not  discussed  but  may  be  inferred  by  the  test  conditions.  Hisiung  and  Chou71  report  an 
increase  in  fractal  dimension  with  increase  in  impact  toughness  in  high  strength  low  alloy 
(HSLA)  specimens  that  should  be  fracturing  through  cleavage  (test  temperatures  ranging 
from  -145  to  -20°C).  Hui  et  al.7-  report  a decrease  in  fractal  dimension  with  increasing 
impact  toughness  in  an  experimental  steel  that  should  be  fracturing  through  microvoid 
coalescence  (test  temperatures  at  -20°C). 

Not  all  studies  using  impact  toughness  are  consistent,  however.  According  to 
Wiencek  and  co-workers,  ' no  relationship  exists  between  the  fractal  dimension  and 
impact  toughness.  In  this  case  it  should  be  mentioned  that  the  authors  used  vertical 
profiles  rather  than  slit-islands.  The  writers  note  that  they  applied  a segment  counting 
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method  similar  to  the  box  counting  method  (Kolmogorov  fractal  dimension).  Despite 
impact  testing  spherodized  steel  specimens  at  temperatures  of-196°C,  -45°C  and  -20°C, 
all  D*  values  were  found  to  be  between  0.09  and  0. 10. 

Relationship  between  fracture  toughness  and  fractal  dimension 

Mu  and  Lung74  have  studied  the  change  in  fractal  dimension  with  fracture 
toughness  in  two  medium  carbon  steels.  In  one  set  of  specimens  the  heat  treatment  was 
altered  to  vary  the  fracture  toughness  while  in  the  second  set  the  temperature  was 
regulated.  In  both  cases  the  study  revealed  a decrease  in  fractal  dimension  with 
increasing  fracture  toughness  as  illustrated  in  Figure  2-10. 


Figure  2-10.  Plot  of  fracture  toughness  against  fractal  dimension,  (a)  The  specimens 
have  been  heat  treated  differently,  (b)  The  temperature  has  been  varied.  [Ref  74]. 

The  basic  relationship  derived  from  this  study  can  be  stated  as 


(2-9)  Ln(Kic)  = Constant+[(  l-DF)Ln(s,)]/2 
Where  the  "constant"  term  incorporates  the  quantities:  E,  y and  v. 
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The  terms  are  defined  as  follows 

Kic=  Plane  strain  fracture  toughness.  E'  = Young's  modulus/!  1-v2). 

Y — Effective  surface  energy  (Ytrue  surface  energy  + YPIastic  strain  energy)- 

Dp  = Surface  fractal  dimension.  v = Poisson's  ratio. 

Sj  = a "step"  of  crack  propagation. 

Mu  and  Lung  prompted  additional  studies  and  many  attempts  to  explain  their 
data.  Many  of  the  studies  employed  impact  toughness  specimens  (see  previous  section) 
which  provide  a measure  of  toughness  but  no  directly  comparable  data.  There  is  no 
confirmed  correlation  between  impact  toughness  and  fracture  toughness.73 

In  follow-on  work  Lung  and  Zhang  explained  that  the  negative  correlation 
between  D*  and  Kic  could  be  explained  in  terms  of  the  intricacies  of  the  perimeter-area 
slit-island  measurement  method  as  well  as  the  micromechanism  of  fracture.76  A 
quantitative  "fractal"  description  of  the  mechanisms  of  transgranular  and  intergranular 
fracture  was  found  to  be  consistent  with  their  qualitative  observations.77 

Most  recently  the  fracture  toughness-fractal  dimension  has  been  examined  by  Su 
and  Lei.  They  attempt  to  explain  all  fracture  toughness-fracture  surface  relationships 
with  a unified  model  based  on  profile  roughness.7*  The  final  form  of  their  equation  is 
(2-10)  J1C  = (aysL)(A0  + B0  ln(AD//) 

Where: 

Jic  = Critical  J-integral  value  at  fracture  initiation. 

oys  = Yield  strength.  L = a microstructural  length.  Ao  and  B()  = Material 

constants 

AD  = D-l  for  profiles,  D-2  for  surfaces. 

/ = volume  fraction  of  dimple  nuclei. 

The  interpretation  in  this  case  being  that  fracture  toughness  depends  on  several 
microstructure  - surface  "roughness"  features  and  that  a single  parameter,  namely  D,  is 
insufficient  to  determine  fracture  toughness  from  fracture  surface  topography  data. 
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Interpretation  of  fractal  dimension  relationships  in  metals 

The  interpretation  of  toughness-fractal  dimension  graphs  derived  from  fracture 
surfaces  have  been  the  focus  of  extensive  research.  This  area  of  study  appears  to  have 
attracted  more  attention  than  fracture  toughness-fracture  surface  roughness  correlations. 
This  may  be  due  to  the  fact  that  researchers  are  still  attempting  to  understand  the  meaning 
of  fractal  curves  and  their  subtle  distinguishing  characteristics. 

One  of  the  significant  observations  reported  by  many  researchers  is  the  reverse 
sigmoidal  shape  (RSC)  of  many  Richardson  or  perimeter-area  type  slit-island 
measurements.  A classic  work  in  this  area  was  published  by  Underwood  and  Banerji.79 
Their  paper  focused  on  the  assertion  that  fractal  plots  were  naturally  reverse  S shaped  and 
that  this  shape  contains  useful  information  about  the  fracture  process.  Underwood  and 
Banerji  approached  this  problem  by  converting  the  RSC  shapes  into  linear  plots  which 
appeared  to  offer  better  resolution  to  differentiate  different  fracture  mechanisms.  The 
technique  was  applied  to  specimens  of  AIS1  4340  steel  heat  treated  and  tempered  at 
different  temperatures.  These  researchers  found  that  the  fractal  dimension  was  lowest 
within  the  tempering  range  which  usually  produces  embrittlement  in  this  type  of  steel. 
Their  key  findings  are  shown  in  Figures  2-1  1 and  2-12.  No  fractographic  information 
discussing  the  fracture  mechanism(s)  observed  was  presented  for  any  of  the  specimens  in 
this  study. 

Other  researchers  have  also  devoted  extensive  effort  into  the  examination  of  this 
often  seen  RSC  behavior.  Some  have  explained  this  fractal  graph  shape  in  terms  of  the 
interaction  of  multiple  fractals  (also  known  as  "multirange  fractals"  which  expose 
different  fractal  behavior  at  different  scales).80  Investigations  by  Shi  et  al.81  suggest  that 
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the  shape  of  fractal  plots  are  influenced  not  only  by  D*,  but  also  by  the  configuration  of 


the  "initiator  length".  As  a result  both  factors  should  be  incorporated  into  any  models 


attempting  to  relate  fracture  toughness  and  fractal  dimension. 
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Figure  2-11.  Fractal  plots  of  AISI  4340 
steel  samples  tempered  at  different 
temperatures  [Ref.  79], 


1.0  5.0  10  50  100  500  1000 

Measuring  unit  (tj), 


Figure  2-12.  (a)  Linearized  fractal  plots,  (b)  Plot  of  fractal  dimension  v. 
tempering  temperature.  According  to  Underwood  and  Banerji  the  lowest  values 
of  D are  found  within  the  temper  embrittlement  region  [Ref.  79], 


RSC  behavior  continue  to  be  a much  disputed  characteristic  of  fractal  graphs. 
Many  researchers  maintain  that  the  inflections  points  have  some  inherent  meaning  and 
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must  be  considered'  ~ while  others  claim  that  the  points  are  an  artifact  of  the  measurement 
method.  Russ4  and  others  have  shown  that  it  is  possible  to  artificially  induce  inflection 
points  at  both  ends  of  the  curve  shown  in  Figure  2-4.  At  the  small  ruler  length  side  of  the 
graph,  the  curve  may  bend  as  a result  of  the  measuring  unit  being  as  small  as  or  smaller 
than  the  profile  width.  At  the  large  ruler  length  side  of  the  plot  the  curve  may  deviate 
from  a straight  line  because  the  ruler  is  approaching  or  surpassed  the  size  of  the  measured 
object  and  the  "step  count"  is  no  longer  accurate.  This  last  explanation  is  the  most 
reasonable.  However,  there  are  conflicting  data. 

It  is  the  view  of  Dauskardt  et  al.83  that  the  former  is  the  case.  That  is,  that  fractal 
curve  inflection  points  reveal  something  intrinsic  about  microstructures,  particularly 
those  at  the  low  ruler  length  side  of  the  plots.  This  team  considered  many  types  of 
metallic  fracture  surfaces  in  a variety  of  steels.  Through  an  analysis  of  artificially 
generated  profiles  and  real  fracture  surfaces  they  conclude  that  fractal  plots  can  indeed 
identify  the  influences  of  microstructure  on  fracture  surface  profiles.  Two  of  their  graphs 
which  may  be  applicable  to  this  study  are  illustrated  in  Figure  2-13. 
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Figure  2-13.  Log-log  R|_  v.  r\  graphs  (a)  and  (b)  illustrate  the  possible  interpretation  of 
fractal  plot  curves  and  their  relationship  to  microstructures.  Note:  RL— Ror)'  'Dl  [Ref  83]. 
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Dauskardt  et  al.  postulate  that  the  slope  fractal  of  graphs  constructed  with  short 
steps  may  be  attributed  to  cleavage  steps  in  brittle  materials  and  slip  steps  in  ductile 
materials.  At  larger  step  lengths  the  plots  are  influenced  by  grain  size  in  brittle  materials 
and  particle  spacing  in  ductile  materials.  Williford  has  also  commented  on  the  perceived 
sigmoidal  shape  of  the  step  size-total  size  fractal  plot.84'86  In  a manner  similar  to 
Dauskart  and  co-workers,  Williford  attributes  this  shape  to  the  influence  of  a variety  of 
fracture  mechanisms  (see  Figure  2-14). 
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Figure  2-14.  Plot  of  step  length  v 
, 4 energy  of  fracture.  Interpretations 

of  various  sections  along  the 
i 3 ~ curve  also  given.  [Refs.  84-86], 
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Based  on  his  observations  in  metals  and  other  materials,  Williford  also  provided  a 
unified  qualitative  model  for  fracture.8^  His  interpretation  proposed  that  materials  may 
be  divided  into  two  basic  categories;  brittle  and  ductile.  In  each  of  these  classes  the 
mechanisms  of  fracture  is  said  to  depend  on  the  characteristics  of  the  microfracturing 
mechanism  prior  to  full  scale  unstable  crack  propagation.  If  the  fracture  mechanism  can 
be  viewed  as  the  nucleation  of  small  cracks,  their  joining  and  propagation,  the  resulting 
fracture  surface  would  be  expected  to  be  rather  flat  and  smooth  at  low  toughness.  As  the 
energy  of  fracture  increases  it  is  reasonable  to  expect  that  the  excess  energy  be  dissipated 
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by  creating  more  surface.  This  process  would  result  in  more  tortuous  fracture  surfaces 
and  resulting  greater  measured  fractal  dimensions. 

Williford  provides  a similar  argument  in  the  case  of  ductile  fracture  surfaces.  In 
this  case,  as  the  fracture  toughness  increases,  the  surface  is  "stretched  out".  That  is,  the 
dimples  that  are  created  tend  to  elongate  and  widen  resulting  in  an  overall  smoother 
surface.  The  proposal  is  illustrated  in  graphical  form  in  Figure  2-15. 


Figure  2-15.  Surface  dimension  (Ds)  vs.  surface 
energy  (y).  [Ref  85] 


Difficulty  correlating  data. 

Williford  made  one  of  the  first  attempts  to  collate  the  known  information  about 
fracture  toughness-fractal  surface  relationships  and  establish  a unified  model.  Despite  the 
noble  attempt,  doing  so  with  the  available  information  may  not  have  been  appropriate. 
Each  of  the  data  sets  referenced  in  Figure  2-15  were  collected  and  analyzed  in  different 
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ways  in  their  respective  sources.  The  data  for  ceramics  was  obtained  by  Mackin  et  al. 
using  the  perimeter-area  relationship  on  slit-islands  and  true  K[c  values,  the  data  for  AISI 
4340  steel  was  collected  by  Underwood  and  Banerji88  using  fracture  profiles  and  a 
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linearization  scheme  for  sigmoidal  curves,  the  data  for  300-grade  maraging  steel  was 
obtained  by  Mandelbrot  et  al.'°by  using  impact  toughness  specimens  and  the  perimeter- 
area  relationship,  and,  finally,  the  data  for  titanium  alloys  was  gathered  by  Pande  et  al.62 
on  dynamic  tear  energy  specimens  (a  measure  of  toughness  per  ASTM  E 604)  using  both 
vertical  sections  and  the  perimeter  area  relationship. 

The  difficulties  encountered  in  the  literature  by  Williford  persist  in  the  present. 
Studies  continue  to  emerge  which  contain  data  that  is  difficult  to  correlate  to  the  work  of 
other  authors.  It  is  one  of  the  aims  of  this  work  to  follow  the  methods  prescribed  by  Hill 
et  al.,14  so  that  the  data  may  be  analyzed  by  researchers  in  quantitative  fractography.  As 
it  stands  today,  the  goal  of  establishing  a definite  microstructure-bulk  material  properties- 
fracture  surface  topography  relationship  or  a model  which  adequately  describes  the 
seemingly  similar  observations  made  in  many  materials  remains  elusive. 


CHAPTER  3 

MATERIALS  AND  METHODS 


Material 

Basis  For  Selection 

A steel  meeting  the  compositional  requirements  of  The  American  Institute  of 
Steel  and  Iron  (AISI)  4340  was  used  in  this  study.  The  alloy  was  selected  on  the  basis  of 
its  ready  availability,  well  documented  mechanical  behavior,  relatively  low  cost,  and 
technological  significance.  The  alloy  is  in  widespread  industrial  use  and  has,  for  decades, 
established  the  foundation  for  the  development  of  many  commercial  high  strength  steels. 

AISI  4340  is  particularly  prevalent  in  the  aerospace  industry.  The  alloy  is 
commonly  found  in  landing  gear  pistons  and  cylinders  as  well  as  many  components  of 
the  flight  control  systems  and  airframes.  If  any  of  these  items  should  fracture,  one  of  the 
first  pieces  of  information  an  investigator  might  want  to  know  is  the  fracture  toughness  at 
the  time  an  incident  occurred.  This  information  may,  in  turn,  be  employed  in 
determining  material  or  environmental  conditions  at  the  time  of  failure. 

This  study  was  pursued  with  the  dual  objectives  of  providing  an  investigator  with 
a useful  tool  to  obtain  this  information  as  well  as  being  able  to  answer  deeper  scientific 
questions  about  the  nature  of  ductile  fracture.  The  use  of  AISI  4340  satisfied  both 
practical  and  scientific  requirements.  If  the  study  were  successful,  it  was  anticipated  that 
the  results  could  be  equally  applicable  to  alloys  that  developed  from  AISI  4340. 
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Alloy  Condition  As  Received 

The  raw  material  was  delivered  as  two  rolled  plates  measuring  1 .52  cm  x 35.7  cm 
x 55.9  cm.  Each  plate  had  been  cut  from  the  same  parent  slab.  A visual  examination 
revealed  no  abnormal  pitting,  cracking  or  unusual  levels  of  corrosion. 

The  chemical  composition  of  one  plate  was  measured  through  Glow  Discharge 
Spectroscopy  (Leco  GDS-750A).  The  results  were  consistent  with  the  requirements  of 
AISI  4340  and  are  given  in  Table  3-1 . 


Table  3-1  Chemical  Composition  of  Raw  Material 


Element 

AISI  4340  (WT%) 

Sample  (WT%) 

Carbon 

0.38  - 0.43 

0.382  -0.412 

Manganese 

0.60  - 0.85 

0.78-0.80 

Silicon 

0.15-0.35 

0.278  - 0.288 

Phosphorus 

< 0.025 

0.0082 

Sulfur 

< 0.025 

0.011  -0.012 

Chromium 

0.70  - 0.90 

0.868  -0.879 

Nickel 

1.65-2.00 

1.687-  1.702 

Molybdenum 

0.20  - 0.30 

0.244  - 0.249 

Copper 

<0.35 

0.143-0.145 

Iron 

94.75  - 96.27 

95.4-95.5 

Fracture  Toughness  Specimens 
Initial  Manufacturing  and  Hardening 

The  two  plates  were  used  to  manufacture  24  rectangular  Compact  Tension  (CT) 
type  fracture  toughness  samples  according  to  the  requirements  of  ASTM  E 399. x4  All 
samples  were  configured  to  meet  a standard  L-T  orientation.  This  designation  has  a 
fracture  plane  whose  normal  is  in  the  longitudinal  direction  (L)  of  the  plate  and  an 
expected  direction  of  crack  propagation  coincident  with  the  width  (also  known  as  the 
long-transverse  or  "T")  direction. 


31 


In  an  effort  to  minimize  cracking  and  machining  difficulties,  the  samples  only  had 
the  pin  holes  drilled  and  its  hole  comers  rounded  prior  to  hardening  (Austenitizing  + Oil 
Quenching  only).  Following  hardening,  the  samples  were  electric  discharge  machined 
(EDM)  to  create  chevron  type  (V-shaped)  starter  notches. 

The  initial  heat  treatment  (with  the  exception  of  tempering)  was  performed  in 
accordance  with  AMS  2759/1, 9,1  Heat  Treatment  of  Carbon  and  Low  alloy  Steel  Parts, 
one  of  the  many  industrial  standards  which  provide  guidelines  for  heat  treating  of  AISI 
4340.  Consistent  with  this  document,  the  samples  were  Austenitized  at  816°C  and  held 
for  a minimum  of  25  minutes.  Again  per  AMS  2759/1,  the  specimens  were  quenched  in 
oil  (at  16°C  - 71°C)  and  cooled  to  room  temperature.  As  a result  of  this  heat  treating 
process,  the  samples  produced  a uniform  through  hardness  in  the  range  of  52  - 56  HRC. 
Precracking 

ASTM  E 399  requires  that  fracture  toughness  samples  have  sharp  cracks  of  a 
prescribed  length  and  shape.  In  this  study,  cracks  of  the  type  required  were  produced  by 
stress  corrosion  cracking  (SCC).  A solution  of  5%  HCl/water  was  applied  by  eyedropper 
to  the  starter  notch  of  a sample.  Five  drops  were  sufficient  to  penetrate  and  wet  the 
notch.  The  area  was  sealed  by  a layer  of  teflon  tape  and  the  specimen  was  loaded  to  2.2 
kN.  Visual  monitoring  revealed  that  this  method  could  produce  starter  cracks  of  the 
desired  length  within  30  minutes.  All  samples  were  pre-cracked  in  the  same  manner. 

Immediately  following  this  operation,  specimens  were  cleaned  and  neutralized  by 
immersion  in  a ultrasonic  cleaning  chamber  filled  with  water/soap  solution  for  30 
minutes.  Specimens  were  then  hot  air  dried  (~  49°C)  for  another  30  minutes  at  50°C  and 


coated  with  a light  oil. 
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Influence  of  SCC  on  Specimen  Fracture  Toughness 

In  high  strength  steels  (hardness  > 40  HRC)  the  mechanism  of  stress  corrosion 
cracking  is  thought  to  act  in  the  same  manner  as  hydrogen  embrittlement.  A fracture  path 
develops  between  prior  austenite  grain  boundaries.  The  key  distinction  between  the  two 
being  the  source  of  hydrogen.  In  SCC  the  hydrogen  is  generated  by  a corrosion  reaction 
while  in  hydrogen  embrittlement  the  source  is  usually  some  processing  operation  such  as 
chrome  or  cadmium  plating. 

The  influence  of  hydrogen  on  plane-strain  fracture  toughness  depends  principally 
on  the  material,  heat  treatment,  quantity  of  hydrogen  present,  and  loading  (including  but 
not  limited  to  rate  and  mode).  All  specimens  in  this  study  were  neutralized,  washed  and 
dried  immediately  after  exposure  to  HC1.  In  addition,  specimens  of  sets  2 through  8 were 
exposed  to  a elevated  tempering  temperature.  Taking  these  two  factors  into  account,  and 
the  fact  that  the  samples  were  allowed  to  rest  for  several  weeks  prior  to  fracture 
toughness  testing,  hydrogen  is  thought  to  have  been  driven  out  of  the  specimens.  The 
contribution  of  hydrogen  to  fracture  toughness  results  is  considered  to  have  been 
negligible  if  any. 

Stress  corrosion  cracking  produces  a different  type  of  crack  tip  geometry  from 
that  obtained  through  ASTM  E399.  The  standard  requires  a fatigue  induced  pre  crack 
produced  under  strict  loading  conditions.  These  condition  often  result  in  a sharp 
transgranular  crack.  SCC  in  these  specimens,  on  the  other  hand,  resulted  in  transgranular 
cracks  also  of  sharp  crack  tip  geometry.  The  fact  that  the  crack  path  and  tips  might  be 
different  between  the  standard  and  the  specimens  was  not  considered  to  be  highly 
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significant.  The  application  of  fracture  toughness  principles  to  cracked  structures,  for 
example,  does  not  require  prior  knowledge  of  fracture  mode. 

Other  reasons  exist  for  considering  SCC  induced  crack  tip  geometry  to  be  neglible 
on  Kic.  Since  these  crack  fronts  are  relatively  straight,  every  downward  tilt  of  a grain 
boundary  must  be  accompanied  by  a upward  grain  boundary  tilt  elsewhere.  Even  though 
a SCC  crack  follows  a more  tortuous  grain  boundary  path,  the  up  and  down  orientations 
of  the  crack  tips  cause  its  effects  to  average  out  to  near  zero  over  the  length  of  a crack 
front. 

Tempering 

The  pre-cracked  collection  of  24  CT  samples  was  divided  into  eight  sets  of  three. 
The  first  set  was  placed  aside  to  form  an  untempered  subset.  The  remainder  were  each 
tempered  at  480°C  in  a industrial  air  furnace.  Each  sample  was  treated  as  indicated  in 
Table  3-2  to  obtain  hardness  values  in  the  range  of  52  to  40  HRC. 


Table  3-2 

Tempering  Schedule  for  CT  Samples 


Set 

Hardness  Before 

Tempering  Treatment 

Hardness  After 

# 

(HRC  ±2) 

Time  & Temperature 

(HRC  ±2.0) 

1 

54 

None 

54 

2 

54 

480°C  for  05  minutes 

52 

3 

54 

480°C  for  10  minutes 

50 

4 

54 

480°C  for  15  minutes 

47 

5 

54 

480°C  for  30  minutes 

46 

6 

54 

480°C  for  50  minutes 

44 

7 

54 

480°C  for  75  minutes 

42 

8 

54 

480°C  for  95  minutes 

40 

Notes:  a)  Samples  were  given  10  minutes  to  come  up  to  temperature  then  soaked  for  the  specified  time. 

b)  Samples  were  immediately  water  quenched  after  tempering. 

c)  Individual  hardness  values  are  given  in  Table  4-1 . 
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Samples  of  40  HRC  were  specifically  desired  due  to  the  technological 
significance  ot  AISI  4340  in  this  form.  Many  aviation  and  industrial  components  are 
used  in  this  hardness  condition  (corresponding  to  a UTS  of  1,241  to  1,379  MPa). 

The  surfaces  of  several  specimens  tempered  for  the  longest  time  were  examined 
under  a SEM  to  determine  the  possible  extent  of  elevated  temperature  oxidation  within 
the  pre  crack  surfaces.  The  size  of  the  oxide  layer  was  estimated  to  be,  on  average,  1 
micron  or  less.  The  thickness  of  the  oxide  layer  is  considered  to  have  had  a negligible 
effect  on  any  subsequent  fracture  toughness  testing. 

Changes  That  Occur  on  Quenching  and  Tempering  AISI  4340  Steel 

A brief  review  of  the  microstructural  changes  that  occur  on  quenching  and 
tempering  AISI  4340  is  warranted  at  this  point.  The  features  that  develop  as  a result  of 
these  processes  significantly  influence  the  overall  mechanical  properties  as  well  as  the 
fracture  process.91 

The  objective  of  most  heat  treatments  involving  medium  carbon-low  alloy  steels 
is  to  produce  a microstructure  consisting  largely  of  tempered  martensite.  Martensite  is 
the  particular  name  assigned  to  the  body  centered  tetragonal  (bet)  structure  which  results 
from  the  diffusionless  transformation  of  face-centered  cubic  (fee)  gamma  iron  (a.k.a. 
austenite)  on  quenching  from  any  temperature  over  727°C.92 

The  usual  austenitizing  temperature  for  AISI  4340  is  approximately  8 1 6°C.  At 
this  temperature,  and  in  the  austenitic  phase,  carbon  may  be  located  along  the  edges  and 
cube  center  of  the  fee  crystal  structure.  Upon  quenching,  the  structure  desires  to 
transform  itself  to  the  equilibrium  bcc  crystal  structure,  however,  the  process  is  inhibited 
from  doing  so  by  the  presence  of  carbon  atoms.  The  resulting  intermediate  bet  unit  cell 
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structure  has  a total  of  four  carbon  atoms  positioned  between  the  iron  atoms  in  the  "c"  or 
vertical  direction  and  an  additional  carbon  atom  at  the  lower  and  upper  square  faces. 

This  carbon-strained  unit  cell  produces  a steel  which  is  very  hard  and  brittle.  At 
this  stage,  the  nearly  100%  martensitic  steel  is  virtually  useless  unless  it  is  softened  and 
toughened  by  tempering.  The  tempering  temperature  and  time  controls  the 
decomposition  of  martensite  as  well  as  the  formation  of  various  other  second  phases. 

In  AISI  4340  steel,  martensite  forms  in  laths  on  quenching.  The  final 
decomposed  microstructure  (generally  achieved  by  long  term  exposure  to  elevated 
temperatures)  is  a mixture  of  iron  and  iron  carbide  in  the  form  of  layers  (pearlitic)  or 
carbide  spheres  in  a iron  matrix  (spherodized).92  In  between  these  two  extremes  there 
may  exist  a multitude  of  intermediate  and  stable  carbides.  These  also  tend  to  change  as  a 
result  of  tempering.  The  driving  forces  for  these  changes  include  a reduction  in  strain 
energy  when  moving  from  bet  to  bcc  structures,  supersaturation  of  carbon  in  martensite,  a 
reduction  in  interfacial  energy  in  going  from  lath  martensite  to  ferrite/cementite 
structures,  and  lowering  of  interfacial  energy  when  coarsening  ferrite/cementite 

92 

structures. 

In  practice,  quenched  steels  having  approximately  0.4%  carbon  also  contain  a 
percentage  of  retained  austenite.  A level  of  austenite  of  about  2 to  4%  by  volume  is  not 
uncommon.  The  decomposition  sequence  of  the  initial  quenched  structure 
(martensite/retained  austenite)  as  it  proceeds  towards  the  equilibrium  structures  (iron/iron 
carbide)  involves  several  stages.  They  are  presently  considered  to  be  as  noted  in  Table  3- 
3. 93 
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The  tempering  steps  reported  in  Table  3-3  are  not  discrete  sequential  steps. 
Considerable  overlap  sometimes  exists  and  it  is  difficult  to  say  specifically  when  one  step 
ends  and  another  begins. 93 


Table  3-3 


Step 

Change 

Controlling  Mechanism 

1 

Retained  austenite  transforms  to  martensite.. 

Diffusionless  transformation. 

2 

Carbon  atoms  redistribute  to  (a)  lattice  defects 
and  (b)  carbon  clusters. 

Volume  diffusion  of  carbon 
atoms. 

3 

Precipitation  of  s and  r|  transition  carbides. 

Diffusion  of  iron  atoms  along 
dislocations. 

4 

Precipitation  of  Haggs  carbide. 

- 

5 

Decomposition  of  remaining  retained  austenite 
to  ferrite/cementite. 

Volume  diffusion  of  carbon 
atoms  in  austenite 

6 

Conversion  of  the  segregated  carbon  and  the 
transition  carbides  to  cementite 

Volume  and  pipe  diffusion  of 
iron  atoms. 

Note:  Summarized  from  ref  93,  pg  669. 


The  evolution  of  the  tempering  process  and  the  particular  steps  completed  at  any 
given  time  depends  on  the  alloy  and  the  temperatures  employed.  The  length  of  time 
required  to  complete  a step  or  the  location  of  an  alloy  along  the  sequence  of  steps  is 
difficult  to  establish  without  the  use  of  sophisticated  equipment  and  time  consuming 
procedures.  Typically,  a thorough  evaluation  would  require  the  combined  use  of 
metallography,  scanning  electron  microscopy,  transmission  electron  microscopy  and 
energy  dispersive  X-ray  spectroscopy. 

A specific  correlation  between  the  tempers  employed  in  this  study  and  their 
corresponding  detailed  microscrotructural  make-up  was  beyond  the  scope  of  this  report. 
This  work  was  necessarily  limited  to  characterizing  some  of  the  basic  features  of  the 
alloy  and  microstructure  employed  in  this  research.  As  a result,  the  reader  is  merely 
advised  that  a host  of  microstructural  changes  which  could  not  be  completely  quantified 
are  likely  to  have  taken  place  during  tempering  and  that  these  changes  undoubtedly 
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influenced  the  fracture  process.  The  specifics  of  how  this  might  have  occurred  is 
considered  throughout  the  Results  and  Discussion  section  of  this  work. 

Second  Phases  in  AISI  4340 

Second  phases  are  among  the  most  important  features  that  develop  on  tempering 
steel.  As  noted  in  Table  3-3,  the  epsilon  (s)  and  p(eta)  transition  carbides  begin  to  form 
early  in  the  tempering  process.  These  carbides  are  difficult  to  resolve  under  a TEM  and 
are  often  discussed  together  as  the  (s/r|)  carbides.  The  crystal  structure  of  epsilon  is 
hexagonal  while  that  of  eta  is  orthorhombic. 

Epsilon  and  eta  carbides  nucleate  on  martensitic  subgrain  boundaries  and/or 
martensitic  lath  boundaries,  grow  by  the  diffusion  of  carbon  and  eventually  transition  to 
cementite.  Prior  to  becoming  cementite  and  at  the  early  stages  of  cementite  growth  these 
particles  are  sub-micron  sized.  The  martensitic  subgrains  are  between  1 .0  to  0. 1 microns, 
or  less,  in  diameter  and  the  grain  boundary  thickness  is  in  the  order  of  20nm.44 

A study  by  Speich44  shows  that  for  a carbon  steel  of  -0.4%  carbon,  tempered 
from  a hardness  of  62  HRC  (760  DPH)  to  a hardness  of  40  HRC  (400  DPEI),  the 
microstructure  has  transitioned  through  steps  1-5  of  Table  3-3.  At  ~40  HRC  the  steel  is 
roughly  at  the  beginning  of  the  ferrite  recovery  and  cementite  spheroidization  process. 
Most  of  the  cementite  nucleated  and  grown  up  to  now  is  in  the  form  of  rods  and  located 
at  martensitic  lath  boundaries. 

AISI  4340  steel  contains  as  its  primary  alloying  elements  nickel,  chrome  and 
molybdenum.  The  central  function  of  these  ingredients  is  to  improve  hardenability, 
however,  chrome  and  molybdenum  are  also  strong  carbide  formers  and  may  also 
participate  in  the  formation  of  second  phases.  Nickel  provides  some  solution 
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strengthening  but  does  not  undergo  any  reactions.  When  tempered  at  below  ~540°C, 
chrome  and  molybdenum  may  enter  the  cementite  phase  by  forming  (Fe,  M)3C  type 
carbides.93 

In  addition  to  the  desirable  martensitic  and  iron  carbide  phases,  AIS1  4340  steel 
may  contain  a wide  variety  of  inclusions.  The  presence,  quantity  and  size  distribution  of 
inclusions  depend  greatly  on  the  steel  manufacturing  method.  In  modem  day  practice, 
foundrys  go  to  great  lengths  to  minimize  inclusion  content,  however,  no  steel  can  be 
made  to  be  completely  clean  in  an  economically  feasible  manner. 

The  most  important  types  of  inclusions  are  those  consisting  of  sulfides,  oxides 
nitrides  and  silicates.9'’  These  basic  compositions  are  MnS,  A1203,  MnAl204,  FeAl204, 
MgAl204,  TiN,  Si02.%  Oxides  are  introduced  into  the  molten  steel  by  adsorption  from 
the  atmosphere,  sulfides  are  present  in  the  initial  iron  ore  and  A1  and  Ti  are  typically 
introduced  as  part  of  the  deoxidation  process. 

Perhaps  the  most  abundant  inclusion  type  is  manganese  sulfide.  Sulfur  is 
naturally  present  in  the  initially  mined  iron  mix  and  is  difficult  to  remove.  The  usual 
manner  of  controlling  the  location  of  this  potentially  harmful  element  in  steel  is  by 
adding  manganese  and  forming  the  far  more  benign  MnS.  This  inclusion  may  solidify  as 
a initially  spheroidal  structure,  however,  hot  working  causes  it  to  assume  a geometry  and 
orientation  consistent  with  the  deformation  process  and  direction. 

Microstructural  Evaluation  Through  Metallography: 

Two  samples  of  the  heat  treated  steel,  one  at  54  HRC  (untempered)  and  another  at 
42  HRC  (tempered.  Set  7 of  Table  3-2),  were  prepared  in  a polymeric  mount  to  reveal 
their  microstructural  features.  The  plane  made  ready  for  examination  was  equivalent  to 
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the  plane  of  crack  propagation  (L-T,  according  to  ASTM  E 399).  The  samples  were 
initially  polished  through  a progression  of  180,  240,  360,  400,  600,  and  800  SiC  grits  in  a 
standard  metallographic  polishing  wheel.  The  wheel  surface  was  changed  to  a cloth  and 
used  to  continue  polishing  through  a progression  of  1 .0  micron  and  0.3  micron  aluminum 
oxide  slurries.  The  last  polishing  step  employed  a cloth  pad  and  a 0.05  micron  diamond 
abrasive  suspension. 

The  polishing  procedure  resulted  in  samples  having  an  essentially  flat,  scratch 
free,  mirror  surface  finish.  The  polished  surfaces  were  then  etched  for  10-15  seconds 
with  a traditional  2%  nitric  acid/water  (a.k.a.  2%  Nital)  etch.  The  etchant  revealed  the 
microstructures  provided  in  Figures  3-1  and  3-2. 

The  microstructural  features  observed  at  both  extremes  of  the  hardness  range 
were  similar.  The  samples  displayed  the  martensitic  laths  typical  of  AISI  4340  steel  heat 
treated  to  40  HRC  and  above.  No  abnormal  structures,  inclusions  or  chemical  anomalies 
(segregation)  were  observed. 


Figure  3-1 . Microstructure  of  a 54  HRC  sample  at  (a)  -400X  and  (b)  -1000X.  Note 
untempered  martensite. 
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Figure  3-2.  Microstructure  of  a 42  HRC  sample  at  (a)  -400X  and  (b)  -1000X.  Note 
tempered  martensite. 

Heat  Treatment  (Quench)  Uniformity 


In  an  effort  to  insure  that  all  samples  exhibited  uniform  heat  treatment  through 
their  thickness,  a microhardness  traverse  was  performed  in  accordance  with  ASTM 

97 

E384,  Knoop  500  g scale.  The  54  HRC  mounted  and  polished  cross  section  employed 
for  microstructural  evaluation  (untempered)  was  used  in  the  traverse.  The  results  of  this 
evaluation  are  shown  in  Figure  3-3. 
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Figure  3-3.  Microhardness  v distance. 

The  mean  microhardness  value  was  measured  to  be  598  HKsoog-  Figure  3-3 
shows  a fluctuating  set  of  points  about  this  mean  value.  The  fluctuations  are  considered 
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normal  for  microhardness  tests.  The  precision  of  the  measurement  is  typically  in  the 
order  of  ±15  HKsoog.  The  initial  quench  is  considered  to  have  produced  a relatively 
uniform  microstructure  and  a homogeneous  set  of  mechanical  properties  through  the 
sample  thickness. 

The  microhardness  test  results  were  found  to  be  consistent  with  the  known  heat 
treatment  response  of  AISI  4340  steel.  Jominy  end  quench  tests  show  steady  hardness 
values  (less  than  1.0  point  variability  on  the  HRC  scale)  for  the  first  19.0  mm  from  the 
quenched  end.91  Considering  that  the  distance  between  the  CT  sample  sides  and  the 
center  is  no  more  than  6.35  mm,  even  hardening  through  the  thickness  should  be 
expected. 

Grain  Size 

Grain  size  measurement  in  hardened  steels  generally  refers  to  the  prior  austenite 
grain  size.  This  measurement  is  notoriously  difficult  to  make  in  these  alloys  due  to  the 
problem  of  differentiating  the  etched  grain  boundaries  from  the  martensitic  laths.  Note, 
for  example,  the  inability  to  distinguish  grain  boundaries  from  the  bulk  microstructure  in 
Figures  3-1  and  3-2. 

ASTM  El  12  is  the  industry  standard  for  grain  size  measurement.  This 
document  suggests  using  a polished  and  etched  metallographic  cross  section  to  obtain 
grain  boundaries  suitable  for  measurement.  Krauss9-  provides  a number  of  methods  for 
revealing  grain  boundaries,  however,  the  techniques  are  cumbersome  and  time 
consuming.  An  alternative  method  is  provided  by  Raymond.99 

Raymond's  method  can  be  performed  on  specimens  where  fracture  surfaces  have 
been  formed  by  intergranular  fracture.  In  this  case,  a set  of  images  may  be  obtained 
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through  a scanning  electron  microscope  (SEM)  from  the  surface.  A series  of  grain 
diameter  measurements  is  collected  and  a chart,  published  in  Raymond's  work,  can  be 
employed  to  obtain  an  equivalent  ASTM  grain  size  measurement. 

Recall  that  stress  corrosion  cracking  (SCC)  was  used  to  create  pre  cracks  in  the 
fracture  toughness  samples.  In  AIS1  4340  SCC  occurs  through  the  prior  austenite  grain 
boundaries.  The  resulting  fracture  surface  is  intergranular  in  nature.  Figures  3-4  (a) 
through  (d)  were  taken  from  the  pre  cracked  areas  of  fracture  toughness  samples  heat 
treated  to  54  HRC  and  left  untempered.  A set  of  sixty  diameter  measurements  were 
taken.  The  mean  diameter  was  measured  to  be  7.7  microns.  This  value  results  in  an 
average  ASTM  grain  size  of  1 1 . 

Microchemistry 

The  surfaces  illustrated  in  Figures  3-1,  3-2  and  3-4  (a)  through  (d)  were  inspected 
under  an  SEM  through  energy-dispersive  X-ray  spectroscopy  (EDS)  in  accordance  with 
the  methods  recommended  in  ASTM  El 508. 100  The  topographies  were  examined  for 
evidence  of  abnormal  composition,  inclusions,  unusual  segregation  or  other  similar 
chemical  anomalies.  The  lenticular  fissures  illustrated  in  Figures  3-4  (b)  and  (d)  were 
given  particular  attention.  Many  similar  crevices  were  observed  throughout  a large 
number  of  the  pre-cracked  fracture  surfaces. 

Figure  3-5  shows  a typical  EDS  spectrum  taken  from  the  intergranular  region  of  a 
fracture  toughness  specimen.  The  spectrum  is  consistent  with  the  composition  reported 
in  Table  3-1 . The  peaks  associated  with  iron,  nickel,  chrome  and  molybdenum  are 
marked.  Elements  present  in  low  concentrations  reveal  peaks  that  are  barely  discernible 
from  the  background  radiation.  The  smaller  peaks  are  expanded  vertically  in  Figure  3-6. 
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Figure  3-4.  Pre  crack  fracture  surfaces  obtained  through  SCC  of  fracture  toughness 
samples  (a  - d).  Granular  features  represent  the  sizes  and  shapes  of  the  prior 
austenite  grains.  The  microstructure  contained  a distribution  of  small  MnS  stringers 
typical  of  AISI  4340  (d)  and  a few  large  stringers  (b). 


Full  Scale  22280  cts  Cursor:  6.390  keV  (24867  cts)  keV 

Figure  3-5.  EDS  Spectrum  of  Overall  Intergranular  Region 
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Figure  3-6.  EDS  Spectrum  of  Overall  Intergranular  Region.  (Spectrum  same  as  Fig. 

3-5  expanded  vertically). 

Figure  3-7  shows  an  EDS  spectrum  typical  of  the  data  collected  from  fracture 
surface  fissures  such  as  those  illustrated  in  Figures  3-4  (b)  and  (d).  Note  the  elevated 
levels  of  manganese  and  sulfur  compared  to  the  overall  spectra  of  Figures  3-5  and  3-6. 
Two  superimposed  and  normalized  spectrums  showing  the  overall  composition  and  that 
found  within  one  of  these  surface  anomalies  is  provided  in  Figure  3-8. 


Figure  3-7.  EDS  Spectrum  from  Fracture  Surface  Fissures. 
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Figure  3-8.  Superimposed  and  Normalized  EDS  Spectrums.  Overall  Fracture  Surface 
(red)  v Fissure  (yellow). 


Tensile  Testing 

Tensile  tests  were  conducted  in  accordance  with  ASTM  E8.  Considering  the  vast 
amount  of  tensile  test  data  available  for  AIS1  4340,  only  three  sets  of  two  samples  each 
were  considered  necessary  to  establish  and  verify  average  properties.  Samples  were 
manufactured  and  heat  treated  to  hardness  values  of  40,  46  and  50  HRC  as  indicated  in 
Table  3-2.  Additional  data  was  obtained  from  a number  of  sources  in  the  literature101'102 
(see  Figure  3-9  for  a summary). 

Fracture  Toughness  Testing 

Fracture  toughness  tests  were  conducted  in  accordance  with  ASTM  E 399.  All 
samples  were  loaded  to  failure  in  an  Instron  Series  IX  Automated  Materials  Testing 
System  equipped  with  a 89  kN  load  cell.  A crosshead  speed  of  0.127  cm/min  was 
selected.  Displacement  was  measured  by  mounting  an  extensometer  at  the  mouth  of  the 
CT  samples.  The  tests  collected  load  displacement  information  as  required  by  the 
standard.  Figure  3-10  shows  representative  CT  fracture  surfaces  after  testing. 
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Hardness  v.  Yield  & Ultimate  Tensile  Strength 


Figure  3-9.  Flardness  v.  yield  and  ultimate  tensile  strengths.  Data 
may  be  used  to  correlate  specimen  hardness  to  average  yield  and 
ultimate  tensile  strengths. 

The  load-displacement  graphs  were  employed  as  indicated  in  ASTM  E 399  to 
determine  the  maximum  load  to  fracture,  Pq  The  fracture  surfaces  were  used  to  measure 
the  size  of  the  crack  (pre-cracked  dimension  "a")  prior  to  fast  fracture.  This  data  can  be 
inserted  into  equation  3-1  to  calculate  a conditional  value  for  Kic  (known  as  Kq).  If  this 
value  was  found  to  be  consistent  with  the  sample  thickness  and  mechanical  properties 
given  by  equation  3-3,  the  test  was  considered  to  be  a valid  plain-strain  fracture 
toughness  test  and  Kq  = Kic.  The  measurements  and  results  discussed  above  are  provided 
in  Table  4-1. 

(3-1)  Kq  = (PQ/BWI/2) /(a/W) 

(3-2)  fla/W)  = (2  + a/W)(0.866  + 4.64a/W  - 13.32a2/W2  + 14.72a3/W3  - 5.6a4/W4) 

(1  -a/Wf2 

(3-3)  B > 2.5  (Kq/cJys)2 
Where 

Kq  = Conditional  stress  intensity  factor.  a = crack  length 

Pq  = Load  to  fracture.  Gys  = Yield  stress 

B = Specimen  thickness 
W - Specimen  width 
/ (a/W)  = Shape  factor  function. 
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Figure  3-10.  Specimen  #21  (a  and  b)  surface  topography  is  typical  of  high  fracture 
toughness  steels.  Specimen  #1  (c  and  d)  surface  topography  is  typical  of  low  fracture 
toughness  steels. 


Fractal  Dimension  Measurements 


Specimen  Preparation 

One  side  of  each  of  the  fracture  surfaces  from  the  fracture  toughness  specimens 
was  used  to  manufacture  a sample  suitable  for  fractal  dimension  measurement.  The  slit- 
island  method  was  selected. 

Preparation  steps  began  by  cutting  the  region  of  each  CT  sample  which  contained 
the  fracture  surface  into  a rectangular  section  measuring  approximately  3.54  cm  x 1 .27 
cm  x 1.5  cm.  Each  section  was  immersed  into  a heated  (88°C)  electroless  nickel  plating 
bath.  A plating  time  of  15  minutes  produced  a uniform  deposit  measuring  6 microns 


thick. 
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Each  of  the  plated  fracture  surfaces  was  carefully  sanded  in  a standard 
metallographic  polishing  wheel  through  a progression  of  400,  600,  800  SiC  grits.  The 
wheel  surface  was  changed  to  a cloth  and  used  to  continue  polishing  through  a 
progression  of  1 .0  micron,  0.3  micron  aluminum  oxide  slurries.  The  last  polishing  step 
employed  a cloth  pad  and  a 0.05  micron  diamond  abrasive  suspension. 

The  polishing  steps  required  extreme  care  and  gentle  applications  of  pressure  to 
avoid  destroying  the  fracture  surface.  This  method  was  successful  in  producing 
numerous  slit-islands  in  all  samples.  The  islands  were  very  small  and  difficult  to  see  by 
the  unaided  eye,  however,  they  could  be  readily  identified  under  an  SEM  (BEI  atomic 
number  contrast  mode). 

An  image  representative  of  a plated,  cut  and  polished  specimen  is  shown  in  Figure 
3-11.  This  sample  is  ready  for  insertion  into  the  SEM. 


Figure  3-11.  CT  specimen  cut,  plated  and 
polished  to  reveal  slit  islands  near  its  center 
(circled). 

Slit-island  Image  Capture 

The  grinding  and  polishing  procedure  described  caused  the  removal  of  nickel 


plate  from  high  spots  on  the  fracture  surfaces.  This  resulted  in  the  production  of 
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numerous  flat  islands  of  steel  surrounded  by  nickel  plate.  These  islands  were  readily 
identified  under  a scanning  electron  microscope  (SEM)  in  backscatter  electron  imaging 
(BEI)  mode. 

A CamScan  MaXim  2040  SL  Scanning  Electron  Microscope  was  applied  to 
capture  the  necessary  images.  The  instrument  was  equipped  with  a LaB6  electron  source. 
The  SEM  was  set  up  to  operate  at  35-40  Kv  and  a working  distance  of  approximately  15 
mm.  The  BEI  detector  employed  was  of  a four  quadrant,  annular,  solid  state  design.  The 
detector  was  positioned  above  the  sample  and  surrounded  the  instrument  pole  piece. 

Three  slit  islands  were  photographed  from  each  sample.  A magnification  of  500- 
800X  was  found  to  produce  sufficient  resolution  and  a clear  view  of  each  area  of  interest. 
The  islands  were  selected  on  the  following  basis 

1.  The  island  must  fill  between  70  to  90  %.  of  the  image  window 

2.  The  island  must  have  an  undamaged  contour  (free  from  pits,  voids,  large 
scratches). 

3.  The  island  must  not  exhibit  metallurgical  irregularities  on  its  perimeter  (free  of 
inclusions  or  MnS  stringers). 

4.  The  image  must  have  good  contrast  between  the  base  metal  and  the 
surrounding  Ni  plate 

5.  Slit  islands  were  to  be  located  within  flat  the  plane-strain  region  of  the  sample 
(away  from  the  angular  shear  lips). 

Images  meeting  the  above  requirements  were  initially  photographed  at  a 
resolution  of  640  x 5 12  in  a Tiff  format.  The  pictures  were  then  converted  to  a JPEG 
format  for  ease  of  transfer  to  the  image  processing  software.  All  islands  are  shown  in 
Appendix  C. 

Image  Processing 

Image  processing  was  performed  in  a standard  desktop  computer  operating  a 
software  program  known  as  Image  Pro  (V.  3.0.01 .00).  Image  Pro  included  an  accessory 
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program  known  as  Materials  Pro  (V.  3.1).  The  latter  software  is  capable  of  performing  a 
large  variety  of  materials  related  image  analysis,  optimization  and  measurements.  The 
combination  of  Image  Pro/Materials  Pro  allowed  for  the  introduction  of  an  image,  semi- 
automatic selection  of  the  island  (through  an  area  of  interest)  and  measurement  of  the 
fractal  dimension. 

In  essence.  Materials  Pro  performs  the  same  operations  an  operator  would 
perform  in  a manually  executed  Richardson  plot.  That  is,  once  the  user  defines  the 
perimeter  to  be  measured,  the  software  uses  a large  number  of  "virtual  yardsticks"  to 
obtain  a plot  of  the  log  distance  v log  stride  length  relationship.104  The  software  then 
calculates  and  inserts  a best  fit  line  into  the  data.  The  decimal  part  of  the  slope  of  this 
plot  is  defined  as  D*;  the  fractal  dimensional  increment. 

Preliminary  Testing  of  the  Fractal  Dimension  Measurement  Method 

As  noted  above,  the  fractal  dimension  is  reported  as  D*  where  D*  is  the  mantissa 
of  the  fractal  dimension  D=2.D*.  This  D*  value  assumes  that  a linear  measurement  (such 
as  that  performed  on  island  perimeters),  i.e.,  1-D,  results  in  the  same  fractal  dimension  as 
a 2-D  measurement.  A series  of  standard  test  images  were  collected.  Their  fractal 
dimension  were  determined  in  an  effort  to  verify  that  the  proposed  measuring  method 
would  provide  accurate  measurements.  These  values  obtained  through  the  test  standard 
were  compared  against  known  reported  or  calculated  data.  The  images  covered  the 
anticipated  range  of  results  from  D*  = 0.000  to  D*  = 0.500  as  indicated  in  Table  3-4. 

The  "textbook"  and  "Koch"  test  images  noted  in  Table  3-4  were  obtained  by 
digitally  photographing  pictures  in  a textbook4  (illustrated  in  Appendix  B).  A similar 
method  would  be  required  in  many  instances  to  obtain  suitable  measurements  from  the 
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actual  sample  surfaces.  Manual  contrast  enhancement  was  then  applied  to  the 
photographs  to  create  images  of  sufficient  contrast  for  boundary  detection  in  the  image 
analysis  system.  A second  digital  image  of  the  first  was  actually  used  for  measurement. 

The  contrast  enhancement  consisted  of  manually  darkening  excessively  light 
areas  and  lightening  excessively  dark  areas.  The  objective  in  each  of  these  cases  was  to 
provide  an  object  which  could  be  readily  detected  and  processed  by  the  software  while 
being  careful  to  preserve  its  original  contour. 

The  fractal  dimension  of  sample  images  la  - lc  were  measured  manually.  Eight 
separate  perimeter  measurements  were  taken  from  each  island  using  different  lengths  of 
measuring  units,  i.e.  "yardsticks".  A graph  of  the  log(measuring  unit)  v.  log(perimeter) 
produced  a linear  relationship.  Automatic  processing  of  the  data  points  by  a curve  fitting 
program  known  as  Curve  Expert100  produced  the  reported  fractal  dimension.  The 
correlation  coefficients  for  the  manual  curve  fits  were  calculated  to  be  in  the  range  of 
0.95  - 0.98.  Based  on  the  test  data,  the  Materials  Pro  image  processing  software  was  able 
to  measure  the  fractal  dimension  to  within  ±6%  of  its  true  value.  The  data  were  more 
consistent  than  the  manual  method  which  could  produce  a discrepancy  of  up  to  1 1%. 

Table  3-4 


Test  Images  and  Resulting  D*  Values 


Image 

Description 

Known 

Calculated 

D* 

Materials  Pro 

% Difference 

Black  circle 

0.000 

0.00 

0 

Textbook 

0.17 

0.18 

+6 

Textbook 

0.23 

0.22 

-4 

Sample  Image  la 

0.19 

0.21 

+ 11 

Sample  Image  1 b 

0.18 

0.19 

+6 

Sample  Image  1 c 

0.20 

0.20 

0 

Koch  Triadic 

0.262 

0.25 

-5 

Koch  Quadratic 

0.500 

0.48 

-4 
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In  obtaining  and  processing  digitized  images  every  effort  was  made  to  control 
imaging  parameters  (brightness,  contrast  and  gamma)  in  such  a way  as  to  insure  that  the 
true  edges  of  the  islands  were  detected  by  the  software.  The  effects  of  each  variable  and 
how  the  perimeter  was  subsequently  detected  could  be  observed  in  real  time  on  a 
computer  monitor.  If  the  effects  were  not  carefully  monitored  on-screen,  adjustments 
within  the  imaging  parameters  could  produce  errors  from  known  values  in  the  fractal 
dimension  of  up  to  ± 10%.  Despite  the  manual  and  imaging  parameter  enhancements,  the 
resulting  data  was  found  to  be  comparable  in  accuracy  to  the  manually  calculated  and 
processed  data. 

Sample  Measurements 

A minimum  of  three  BEI  images  were  obtained  from  each  of  the  twenty  four 
samples.  The  images  are  provided  in  Appendix  C.  Each  image  was  used  "as  is"  if  the 
image  analysis  program  could  accurately  distinguish  its  boundaries  from  the  surrounding 
nickel  plate.  Otherwise,  the  island  was  printed  on  paper  for  manual  contrast 
enhancement,  re-photographed  and  digitally  analyzed.  Fractal  dimension  results  are 
provided  in  Table  3-5. 
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Table  3-5 

Fractal  Dimension  Results 


D*  for  Individual  Images 

Sample 

# 

Image  a 

Image  b 

Image  c 

Image  d 

Image  e 

Image  f 

Mean 

D* 

1 

0.20918 

0.18508 

0.20299 

0.20 

2 

0.21019 

0.21092 

0.18999 

0.20 

3 

0.17086 

0.16807 

0.22111 

0.19 

4 

0.31824 

0.28206 

0.24483 

0.28 

5 

0.29884 

0.27498 

0.21906 

0.26 

6 

0.26536 

0.29401 

0.32258 

0.29 

7 

0.24242 

0.22037 

0.30174 

0.23678 

0.25 

8 

0.27602 

0.22414 

0.22633 

0.24019 

0.21861 

0.24 

9 

0.33127 

0.23291 

0.21328 

0.24006 

0.28716 

0.23661 

0.26 

10 

0.21310 

0.19569 

0.17148 

0.19 

1 1 

0.27696 

0.20276 

0.17871 

0.22 

12 

0.21841 

0.21836 

0.19293 

0.21 

13 

0.20662 

0.16948 

0.15626 

0.18 

14 

0.22060 

0.21511 

0.22641 

0.22 

15 

0.21770 

0.15615 

0.18659 

0.19 

16 

0.18952 

0.12210 

0.15969 

0.16 

17 

0.12514 

0.10469 

0.13395 

0.12 

18 

0.11616 

0.14686 

0.13636 

0.13 

19 

0.13191 

0.14458 

0.07253 

0.12 

20 

0.12714 

0.13013 

0.10179 

0.12 

21 

0.12966 

0.08658 

0.07294 

0.11991 

0.04018 

0.09440 

0.09 

22 

0.15203 

0.18989 

0.18587 

0.18 

23 

0.09862 

0.09889 

0.18028 

0.13 

24 

0.13333 

0.14525 

0.19721 

0.16 

Note:  Std  deviation  in  D*  provided  in  Table  4-1 . 


CHAPTER  4 

RESULTS  AND  DISCUSSION 


Tabular  and  Graphical  Summary 

Table  4-1  and  Figure  4-1  summarize  the  results  obtained  through  fracture 
toughness  and  fractal  dimension  measurements. 


Table  4-1  Fracture  Toughness  - Fractal  Dimension  Results. 


Set  & 
Sample 
Numbers 

Hardness 

Bef./Aftr. 

(HRC) 

Crack 

Length 

(m) 

Fracture 

Load 

(Newtons) 

Shape 

Factor 

Kic 

(MPa  Vm) 

D* 

STD  Dev 
in  D* 

#1)  1 

56.3/56.3 

0.02667 

6978.9 

10.45 

25.6 

0.20 

0.010 

2 

56.4/56.4 

0.02540 

8153.2 

9.66 

27.6 

0.20 

0.010 

3 

55.7/55.7 

0.02692 

7116.8 

10.63 

26.5 

0.19 

0.024 

#2)  4 

55.4/52.4 

0.02540 

7561.6 

9.66 

25.6 

0.28 

0.030 

5 

55.5/52.0 

0.02464 

8495.7 

9.23 

27.5 

0.26 

0.033 

6 

55.7/51.7 

0.02464 

9340.8 

9.23 

30.2 

0.29 

0.023 

#3)  7 

53.9/49.7 

0.02438 

16680.0 

9.09 

53.2 

0.25 

0.031 

8 

55.8/49.6 

0.02540 

14233.6 

9.66 

48.2 

0.24 

0.021 

9 

56.1/49.9 

0.02286 

21350.4 

8.34 

62.4 

0.26 

0.040 

#4)  10 

52.5/46.8 

0.02540 

16457.6 

9.66 

55.7 

0.19 

0.017 

11 

51.5/48.0 

0.02616 

16457.6 

10.12 

58.4 

0.22 

0.042 

12 

52.8/47.5 

0.02464 

17792.0 

9.23 

57.6 

0.21 

0.012 

#5)  13 

52.5/45.7 

0.02286 

28467.2 

8.34 

83.3 

0.18 

0.021 

14 

52.3/46.5 

0.02426 

21795.2 

8.98 

68.6 

0.22 

0.005 

15 

51.5/46.8 

0.02426 

21128.0 

8.98 

66.5 

0.19 

0.025 

#6)  16 

51.8/43.0 

0.02464 

30691.2 

9.23 

99.3 

0.16 

0.028 

17 

52.5/43.8 

0.02540 

28467.2 

9.66 

96.4 

0.12 

0.012 

18 

53.1/44.2 

0.02718 

22684.8 

10.80 

85.9 

0.13 

0.013 

#7)  19 

52.6/41.3 

0.02451 

36473.6 

9.09 

116.3 

0.12 

0.031 

20 

52.0/42.0 

0.02489 

36028.8 

9.37 

118.4 

0.12 

0.013 

21 

51.4/42.1 

0.02731 

27724.4 

10.89 

105.9 

0.09 

0.030 

#8)22 

55.5/40.2 

0.02540 

35584.0 

9.66 

120.5 

0.18 

0.017 

23 

50.6/40.5 

0.02540 

35806.4 

9.66 

121.3 

0.13 

0.038 

24 

52.3/40.4 

0.02464 

35673.0 

9.23 

115.5 

0.16 

0.028 
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Fracture  Toughness  v D* 


Figure  4-1 . Fracture  Toughness  v D*. 


The  results  illustrated  in  Figure  4-1  show  a decrease  in  D*  with  increasing 


fracture  toughness.  Some  inconsistencies,  however,  were  noted  in  the  extremes  of  the 


curve.  At  the  low  fracture  toughness  values,  the  D*  values  dropped  unexpectedly.  At  the 


high  fracture  toughness  end,  the  D*  were  observed  to  rise  above  the  downward  trend. 


Fractography 

Specimens  from  all  fracture  toughness  groups  were  selected  and  prepared  for 
fractographic  examination  under  a scanning  electron  microscope  (SEM).  This  inspection 
was  performed  in  an  effort  to  identify  the  primary  fracture  mechanism  for  each  fracture 
toughness  range.  Images  representative  of  the  various  fracture  surfaces  examined  are 
presented  in  Figure  4-2  (a)  through  (f). 
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(d)  Specimen  #13;  83.3  MPa. Vm. 


(f)  Specimen  #23;  121.3  MPa. Vm. 


Figure  4-2.  Images  are  representative  of  the  fractographic  features  observed  in  the 
fracture  toughness  specimens,  (a)  shows  predominantly  intergranular  separation  and 
cleavage,  (b-d)  show  a mixed-mode  fracture  surface  characterized  by  a blend  of 
intergranular  separation,  cleavage  and  microvoid  coalescence,  (e-f)  are  typical  of 
specimens  of  46  HRC  and  below,  are  characterized  by  nearly  100%  microvoid 
coalescence. 
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Figure  4-2  (a)  is  typical  of  untempered  fracture  toughness  specimens  (56  HRC). 
These  surfaces  revealed  predominantly  inter  and  transgranular  fracture  paths  with  little 
ductility.  The  transgranular  appearance  is  also  known  as  cleavage  fracture.105  Some 
limited  ductility  in  the  form  of  dimple  rupture  was  observed  within  grain  boundaries. 
Measurement  results  of  the  percentage  of  surface  area  covered  by  microvoids  in  this  and 
the  remaining  samples  are  reported  in  the  next  section. 

Figures  4-2  (b)  through  (d)  are  typical  of  lightly  tempered  fracture  toughness 
specimens  (52  and  49  FfRC,  respectively).  These  surfaces  revealed  a blend  of 
intergranular  separation,  cleavage  and  dimple  rupture.  When  cleavage  and  dimples  are 
observed  together  this  type  of  fracture  pattern  may  also  be  referred  to  as  quasi-cleavage. 
The  percentage  of  surface  area  covered  by  microvoids  increased  with  increasing 
tempering  time  and  lower  hardness. 

Figures  4-2  (e)  and  (f)  are  typical  of  moderately  to  highly  tempered  fracture 
toughness  specimens  (47  HRC  or  lower)  These  samples  revealed  surfaces  dominated  by 
microvoid  coalescence.  Samples  at  the  lower  end  of  the  hardness  range  revealed  nearly 
100%  dimple  rupture. 

Fracture  Mechanisms 

Scanning  electron  microscopy  of  the  fracture  surfaces  showed  one,  two  or  three 
fracture  mechanisms.  Samples  were  observed  to  have  separated  through  the  grain 
boundaries  (also  known  as  intergranular  or  "IG"  separation),  cleavage  (transgranular) 
and/or  microvoid  coalescence.  These  mechanisms  are  related  to  the  degree  of  ductility 
available  within  the  specimens.  Steel  specimens  of  little  ductility  are  known  to  fracture 
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through  their  grain  boundaries  and/or  cleavage.  In  these  cases,  the  crack  front  essentially 
propagates  through  or  across  the  grains  with  little  to  no  plastic  deformation.10'' 

Microvoid  coalescence  is  typically  the  dominant  fracture  mechanism  in  metallic 
alloys  that  exhibit  moderate  to  high  ductility.  Voids  begin  to  form  in  the  plastic  zone 
ahead  of  a front.  Dimples  in  steel  will  nucleate  preferentially  at  sites  which  contain 
material  defects  or  inhomogeneities  such  as  inclusions,  second-phases,  grain  boundaries 
and  dislocation  pile-ups.105  As  the  strain  in  the  plastic  zone  increases,  the  dimples  grow 
and  eventually  join  to  create  a continuous  fracture  surface. 

Three  images  were  recorded  for  each  fracture  toughness  set  for  measurement  of 
their  respective  proportions  of  intergranular  separation,  cleavage  fracture  and  ductile 
fracture  (a  general  term  which  includes  dimple  rupture,  tearing  and  localized  plastic 
flow).  Once  each  image  was  collected,  it  was  printed,  manually  enhanced  to  readily 
distinguish  ductile  areas  from  cleavage  and  intergranular  facets,  rephotographed  and 
computer  processed. 

Computer  processing  was  accomplished  by  employing  the  thresholding  function 
provided  by  the  "feature  measurements"  subroutine  available  within  Matrials  Pro.  This 
function  is  able  to  inspect  an  acquired  image  and  determine  the  proportion  of  colored 
areas  to  non-colored  areas  as  determined  by  a user.  In  this  manner,  the  percentage  of 
surface  area  covered  by  ductile  fracture  mechanisms  (as  manually  colored,  or  defined,  by 
a operator)  may  be  determined.  The  results  are  provided  in  Table  4-2. 

The  data  presented  in  Table  4-2  was  compared  against  the  fracture  toughness  - 
fractal  dimension  results  of  Figure  4-1 . One  of  the  possible  explanations  for  the  drop  in 
fractal  dimension  values  at  low  fracture  toughness  (untempered  specimens  of  Set  # 1 ) is 
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that  the  fracture  mechanism  in  these  specimens  is  dominated  by  IG  separation  and 
cleavage.  At  the  higher  fracture  toughness  values  (above  -90  MPaVm)  the  possible 
explanation  for  the  deviation  in  linearity  may  be  related  to  the  departure  from  ideal 
linear-elastic  conditions.  The  strict  sample  thickness  requirements  ("B"  dimension)  of 
ASTM  E399,  as  provided  in  equation  3-3,  were  no  longer  satisfied.  At  above  -90 
MPaVm  the  data  should  be  considered  to  represent  a Kq  value  rather  than  a Kic.  The  load 
displacement  diagrams  for  the  high  fracture  toughness,  non-compliant,  E-399  specimens 
is  provided  in  Appendix  A. 


Table  4-2.  Percentage  of  Area  Covered  by  Ductile  Fracture  v.  Sample  Set. 


Sample  Sets  & Mean  HRC 

% Area  Covered  by  Ductile  Fracture 

Set#  1 (56.1  HRC) 

-17% 

Set  # 2 (52.0  HRC) 

-36  % 

Set  # 3 (49.7  HRC) 

-50  % 

Set  # 4 (47.7  HRC) 

-77  % 

Set  # 5 (46.3  HRC) 

100% 

Set  # 6 ( 43.6  HRC) 

100% 

Set#  7 (41.8  HRC) 

100% 

Set  # 8 ( 40.4  HRC) 

100% 

Fracture  Mechanism  v.  Fracture  Toughness  and  D* 


The  fracture  mechanism  observations  made  through  scanning  electron  microscopy 
and  their  relationship  to  the  fracture  toughness  - D*  data  are  presented  in  Figure  4-3. 

This  figure  illustrates  that  the  set  of  24  fracture  toughness  samples  can  be  divided  into 
essentially  two  categories  based  on  the  predominant  fracture  mechanism.  The  first 
category  is  composed  of  samples  where  the  fracture  mode  resulted  in  approximately  83% 
or  more  IG  and  cleavage.  The  second  category  consists  of  samples  where  IG  separation 
and  cleavage  covered  approximately  64%  or  less  of  the  fracture  surface.  Samples  in  the 
second  category  appeared  to  follow  a consistent,  nearly  linear,  pattern. 
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Fracture  Toughness  v D* 


Figure  4-3.  Fracture  Toughness  v.  D*.  Fractographic  features  and  validity  of  fracture 
toughness  results  are  noted. 


Fracture  Surface  Morphology  of  Different  Regions. 

The  slit  islands  used  to  obtain  D*  from  all  CT  samples  were  collected  from  the 


1/3  to  2/3-unstable  crack  regions  of  their  respective  fracture  surfaces  (see  Figure  3-11). 


Though  this  is  not  the  area  where  unstable  fracture  initiates,  its  features  are  thought  to  be 
related  to  the  material  toughness  properties.14  Slit  islands  images  were  recorded 
essentially  at  random  within  these  general  areas.  In  some  specimens  the  islands  were 
obtained  primarily  from  the  1/3  crack-length  region  while  in  others  the  images  were 
obtained  from  the  latter  2/3  crack-length  region.  Observations  were  made  under  a SEM 
to  determine  if  this  selection  process  may  have  influenced  the  character  of  the  data.  In 
other  words,  were  there  any  significant  changes  throughout  the  fracture  surfaces  which 
could  account  for  any  differences  in  D*  results?  No  visible  differences  were  identified 
between  these  regions. 
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Crack  Tip  v.  Bulk  Fracture  Surface  Morphology. 

One  fracture  surface  from  each  of  Sets  1 through  8 were  inspected  at 
magnifications  of  800  to  3,000X.  Examination  of  the  early  regions  and  comparison 
against  the  later  regions  of  the  fracture  surfaces  revealed  few  discernible  differences. 

One  distinction  discovered  to  be  true  for  most  specimens  is  that  the  material  immediately 
adjacent  to  the  crack  tip  was  found  to  be  more  deformed  (stretched)  than  the  bulk  of  the 
fracture  surface.  This  is  typical  of  ductile  surfaces  and  attributed  to  the  uniaxial/biaxial 
stress  state  existing  right  at  the  crack  tip  (strains  are  very  high  and  transverse  stresses  go 
to  zero  at  the  free  surface).  Slightly  inward  from  the  stretched  zone,  yet  still  within  1 to 
2§ic  of  the  crack  tip,  the  surfaces  appeared  to  be  identical  to  the  remainder  of  the  fracture 
surface.  Images  representative  of  the  crack  tips  and  central  sections  typical  of 
intermediate  and  high  fracture  toughness  samples  are  provided  in  Figures  4-4  to  4-7.  For 
these  specimens  crack  tip  opening  displacements  (8ic)  were  approximately  5 pm  (#  10) 
and  1 5 pm  (#  1 7).  Further  discussion  of  the  findings  requires  prior  knowledge  of  ductile 
fracture  in  steels. 

Ductile  Fracture  in  AISI  4340  Steel 

Overview 

Ductile  fracture  is  defined  as  a failure  mechanism  which  takes  place  with 
considerable  plastic  deformation  and  the  consumption  of  large  amounts  of  energy.  In 
most  metallic  alloys,  including  AISI  4340  Steel,  ductile  fracture  is  widely  regarded  as  a 
three  step  process  consisting  of  void  nucleation,  growth  and  coalescence.105  All  three 
processes  together  lead  to  the  formation  of  fracture  surfaces  and  are  influential  in 
establishing  many  alloy  properties  such  as  tensile  ductility  and  fracture  toughness. 
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Figure  4-4.  Crack  tip  region  of  specimen  #10  (Kic  = 55.7  MPa  Vm). 
Transition  region  from  IG  fracture  to  microvoid  coalescence  is  evident. 


Figure  4-5.  Central  fracture  surface  region  of  specimen  #10  (Kic  = 55.7 
MPa  Vm).  ~77  % ductile  fracture. 
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Transition  region  from  IG  fracture  to  microvoid  coalescence  is  evident. 


Figure  4-7.  Central  fracture  surface  region  of  specimen  #17  (Kic  = 96.4 
MPa  Vm).  100  % microvoid  coalescence  is  evident. 


64 


The  variables  which  control  the  void  formation  process  exert  a large  influence  on  ductile 
fracture. 

While  there  are  thee  identifiable  steps  in  this  process,  they  are  not  separable, 
discrete,  events.  Cavities  of  different  sizes  will  nucleate,  grow  and  join  at  different  times 
during  the  deformation  span. 

Void  formation  and  its  influence  on  mechanical  properties  is  an  area  of  intense 
study.  A large  number  of  features  with  resulting  effects  on  properties  have  been 
identified.  Microvoids  are  known  to  nucleate  around  material  discontinuities  such  as 
inclusions,  second  phase  particles,  pores,  grain  or  twin  boundaries  and  dislocation  pile- 
ups.106  Each  of  these  material  characteristics,  when  present,  will  exercise  its  own 
particular  influence  in  the  overall  mechanical  properties. 

Stress  state  is  an  equally  important  factor  in  void  development.  These  features  are 
favored  by  the  presence  of  a high  triaxial  state  of  stress  (high  triaxiality).  Studies  show 
that  under  low  triaxiality  conditions,  void  nucleation  and  growth  may  be  delayed  to 
higher  strains  or  suppressed  entirely. 107  A high  triaxial  stress  state  is  present  in  front  of  a 
highly  loaded  crack  tip  due  to  a thickness  induced  effect  known  as  constraint  (plane- 
strain  conditions).  The  degree  of  triaxiality  if  often  given  as  the  unitless  ratio  of  mean 
(hydrostatic)  to  the  effective  von-Mises  stress  (am/aef).108 
The  Influence  of  Second  Phases. 

In  general,  and  in  the  absence  of  gross  material  flaws  (microcracks  or  pores),  void 
nucleation  begins  first  at  the  largest  inclusions  or  second  phases.  In  steel  tensile 
specimens  this  will  occur  at  relatively  low  strains  slightly  beyond  the  yield  stress.106  The 
larger  the  second  phase,  the  lower  the  strain  required  for  void  initiation. 
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A number  of  void  related  factors  and  features  are  involved  in  determining  when 
nucleation  will  take  place.  Among  these  is  the  condition  of  the  second  phase  (cracked  or 
uncracked),  phase  characteristics  (brittle  or  ductile),  phase  shape  (elongated,  round, 
jagged  etc...),  interfacial  binding  strength  (weak  or  strong),  and  phase  orientation  (long 
axis  oriented  with  or  perpendicular  to  the  loading  axis).106 

In  addition  to  geometry  effects,  the  overall  content  and  distribution  in  space  must 
be  considered.  The  particles  may  be  present  at  a large  or  small  volume  fraction.  Also, 
for  constant  volume  fractions,  the  phase  may  be  large  and  widely  dispersed  or  small  and 
close  together.106 

When  second  phases  are  very  small  (sub-micron  sized),  voids  may  not  nucleate 
from  them  at  all.1 10  In  these  cases,  nucleation  occurs  at  the  large  particles  and  the  smaller 
ones  serve  to  alter  the  flow  properties  of  the  matrix.  In  most  cases,  very  small  particles 
will  provide  a strengthening  effect  consistent  with  the  degree  of  lattice  strain  it  introduces 
or  the  resistance  to  dislocation  motion  that  it  creates. 

Microvoid  Coalescence  in  AISI  4340  Steel. 

Void  formation  in  AISI  4340  steels  heat  treated  to  -1378  MPa  proceeds  first  by 
the  nucleation  and  growth  of  cavities  from  MnS  phases  and  other  similar  large  non- 
metallic  inclusions  (primary  voids).  Once  the  void  size  has  reached  a critical  value,  the 
cavities  link-up  (coalesce)  by  a process  known  as  void  sheet  formation.1 1(1  In  effect,  the 
cementite  particles  in  the  highly  strained  regions  between  large  inclusions  nucleate  and 
grow  their  own  voids  (secondary  voids).  The  combined  effects  of  MnS  nucleated  voids 
connected  by  cementite  nucleated  void  sheets  results  in  a continuous  fracture  surface 
characterized  by  a wide  distribution  of  microvoid  sizes. 
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The  relatively  straight  forward  process  of  microvoid  coalescence  is  controlled  by 
a wide  variety  of  factors.  As  previously  noted,  the  characteristics  of  the  particle/void 
interaction,  as  well  as  the  distances  between  voids,  substantially  affect  the  process.  A 
significant  influence  is  also  exerted  by  the  matrix  properties. 

Fracture  Toughness  of  Steels  Based  on  Microstructural  Features. 

Garrison  et.  al.  have  studied  the  influence  of  second  phases  on  the  void  formation 
process  and  fracture  toughness.1 11 '' 19  This  research  group  began  its  work  with  the  notion 
that  the  Rice  and  Johnson1-0  relationship  between  the  critical  crack  tip  opening 
displacement  (5ic  - a measure  of  toughness)  and  particle  spacing/particle  radius  ratio 
(Xo/Ro)  should  be  modified  to  account  for  alloy  flow  properties.  With  this  in  mind. 
Garrison  proposed  that  5IC  should  scale  with  Xo/CRv/ROk1  1 1 where  Rv  is  a void  radius 
and  R,  is  a inclusion  radius  and  the  evaluation  is  performed  at  an  average  inclusion  radius 
Ro  (Rv/R,  represents  a measure  of  microstructural  ductility).  Their  work  has  found  that  a 
reasonable  correlation  exist  between  primary  void  initiating  particle  spacing  and  the  void 
growth  ratio  as  shown  in  Figure  4-8.  As  the  researchers  indicate,  however,  the 
correlation  is  limited  only  to  alloys  of  "small  particle  spacing  and  limited  ductility".1 15 
At  very  high  toughness  values,  which  vary  according  to  the  individual  alloy,  the  linear 
relationship  no  longer  holds. 

The  notion  that  fracture  toughness  is  related  to  particle  spacing  and  matrix 
properties  has  existed  for  some  time.  An  analysis  by  McClintock121  for  a strain 
hardening  alloy  containing  cylindrical  holes  suggests 

(4-1)  8f  = ln(dp/Dr)(  1-n) 

sinh[(  l-n)(aa  + ob)/(2avm/V3)] 
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where 

Sf*  = Local  fracture  strain. 

dp  = Void  initiating  particle  spacing. 

Dp  = Void  initiating  particle  diameter, 
n = Strain  hardening  exponent. 
oa  = Transverse  stress. 

Ob  = Transverse  stress. 
avm  = Von  Mises  stress. 

Note:  (aa  + ab)/(2avm/V3)  represents  the  local  stress  state. 
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Figure  4-8.  Correlation  between  CTOD  and  particle 
spacing-void  growth  ratios  for  a number  of  steels. 

[ref.  101]. 

While  equation  4-1  is  does  not  conform  with  most  fracture  data,  it  is  considered  to 
be  consistent  with  the  observed  material  trends.  That  is,  if  void  fraction,  strain  hardening 
and  triaxiality  are  made  to  increase  the  local  fracture  strain  will  decrease.  Fracture  strain 
scales  closely  with  fracture  toughness. 

D*  in  High  Fracture  Toughness  Specimens 

Figure  4-3  shows  that  specimens  of  high  fracture  toughness  do  not  appear  to 
follow  the  same  downward  trend  with  increasing  K|C  (or  Kq)  as  others  in  the  same  mixed 
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mode  or  dimple  rupture  fracture  mechanism  category.  That  is,  at  high  fracture  toughness, 
some  samples  appeared  to  increase  in  D*. 

Deviations  from  linearity  may  be  entirely  attributable  to  differences  in  stress  state. 
Note  that  the  highest  fracture  toughness  specimens  did  not  meet  the  sample  dimensional 
requirements  of  ASTM  E399.  The  undersize  samples  (no  longer  under  plane-strain 
conditions)  result  in  a decrease  in  crack  tip  constraint.  Changes  in  constraint  are  known 
to  change  void  growth  rates  and  proportions  (width/depth  as  well  as  Rv/R,).109’ 1 10  When 
triaxiality  is  lowered  void  formation  is  inhibited  and  the  cavities  do  not  grow  as  wide. 

Comparison  Against  Other  D*/Fracture  Toughness  Results. 

The  results  obtained  in  the  present  study  have  been  observed  in  several  previous 
studies.  A significant  problem  in  making  these  comparisons,  however,  is  that  many 
earlier  studies  did  not  measure  fracture  toughness  directly  but  rather  measured  impact 
toughness.  Empirical  relationships  exists  between  these  two  values,  however,  the 
equations  are  highly  dependent  on  the  individual  alloys  and  heat  treatments.  With  the 
exception  of  a fairly  limited  class  of  extensively  studied  alloys  (mostly  ASTM  type 
pressure  vessel  and  bridge  steels),  impact  specimen  values  are  not  generally  used  in  place 
of  actual  fracture  toughness  data.  It  is  accepted,  nevertheless,  that  both  of  these  measures 
of  toughness  will  follow  a similar  trend  as  other  material  properties  are  altered.75' 122 

In  the  case  of  steels,  the  studies  that  are  most  similar  to  this  one  are  those  by 
Mandelbrot,  Passoja  and  Paullay,10  Hilders  and  Pilo,70  Ray  and  Mandal69  and  Mu  and 
Lung.74  The  essence  of  these  studies  have  been  discussed  in  Chapter  2 and  presented  in 
Figures  2-9  and  2-10.  When  compared  to  those  studies,  the  results  presented  here  are 
observed  to  be  consistent.  There  is  an  increase  in  D*  with  an  increase  in  Kic  for  steels  in 
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a brittle  condition  (fracture  primarily  through  cleavage)  and  a decrease  in  D*  with  an 
increase  in  Kic  for  samples  in  a ductile  condition  (fracture  primarily  through  microvoid 
coalescence).  The  results  of  this  study  are  also  consistent  with  the  observations  presented 
by  Williford*5  and  repeated  in  Figure  2-15.  The  data  for  the  AISI  4340  steel  specimens 
agree  with  the  proposal  that  D*  does  not  purely  increase  or  decrease  with  fracture  energy. 
Instead,  it  is  apparent  that  fracture  surfaces,  including  those  of  this  study,  scale  differently 
depending  on  whether  the  materials  behave  in  a primarily  brittle  or  ductile  manner. 
Mathematical  Curve  Fitting  Models. 

An  examination  of  the  data  collected  in  this  study  shows  that  a number  of  linear 
or  nearly  linear  mathematical  descriptions  would  be  appropriate.  The  literature  provides 
few  suggestions  as  to  what  functions  might  serve  to  establish  a suitable  K!C  v D* 
relationship.  Most  studies  with  metals  have  explored  the  fractal  qualities  of  fracture 
surfaces  but  have  failed  to  offer  mathematical  relationships  between  these  quantities  or 
explanations  for  why  this  behavior  might  exist.  The  two  most  promising  equations  are 
provided  by  Mu  and  Lung74  and  Mecholsky  et  al.30  These  functions  are 

4-2)  Ln(Kic)  = Constant+[(l-DF)Ln(si)]/2  (Mu  and  Lung) 
where  the  "constant"  term  incorporates  the  quantities:  E',  y and  v. 

The  terms  are  defined  as  follows 
Constant  = l/2[ln(2E'y)] 

KiC=  Plane  strain  fracture  toughness.  E'  = Young's  modulus/(  1-v2). 

Y — Effective  surface  energy  (ytrue  surface  energy  YPlastic  strain  energy)- 

DF  = Surface  fractal  dimension.  v = Poisson's  ratio. 

8,  =a  "step"  of  crack  propagation  (considered  by  Lung  & Mu  to  be  -constant). 

and 

(4-3)  K|C  = E(a<)D*)l/2 


(Mecholsky  et  al.) 
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where 

Kic  = Plane  strain  fracture  toughness.  E = Young's  modulus, 
ao  = A material  dependent  structure  parameter. 

D = Surface  fractal  dimension. 

Comparison  Against  Lung  and  Mu's  Data: 

The  data  obtained  by  Lung  and  Mu  is  presented  with  the  data  from  this  study  in 
Figure  4-9.  The  format  is  identical  to  that  suggested  by  these  authors  in  their  paper.  The 
data  of  this  study  compares  favorably  and  behaves  in  a similar  manner.  The  appropriate 
functions  and  R2  values  are  noted  on  the  graph. 


Ln(Klc)  v D* 


5 

4 

*75 

* 3 
c 

_i 

2 

1 

0 


Y = -7.59X  + 5.51 
R2  = 0.94 


Y = -4.72X  + 5.55 
R2  = 0.96 


Y = -7.26X  + 5.62 
R2  = 0.79 


♦ Lung  & Mu.  (Tmp  changed) 
■ Lung  & Mu.  (H.T.  Changed) 
▲ This  study  (H.T.  changed) 


0.05  0.1  0.15  0.2 

D* 


0.25 


0.3 


0.35 


Figure  4-9.  Graph  of  Ln(Kic)  v.  D* 


In  applying  the  data  in  this  study,  the  numbers  corresponding  to  the  samples 
exhibiting  cleavage  fracture  were  not  used.  This  set  (set  #1 ) was  excluded  on  the  basis 
that  the  fracture  mechanism  is  significantly  different  from  those  exhibiting  large 
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proportions  of  dimple  rupture.  The  data  belonging  to  set  #1  is  considered  to  behave  more 
like  the  brittle  specimens  of  other  studies. 

Metallic  and  ceramic  samples  which  exhibit  little  ductility  seem  to  have  an 
increasing  value  of  D*  with  increasing  K.[C.  The  three  lowest-toughness  samples  in  this 
study  (Fig.  4-8)  appear  to  follow  this  trend.  As  such,  their  inclusion  within  a larger  set  of 
samples  which  exhibit  primarily  ductile  fracture  mechanisms  that  show  a decrease  in  D* 
with  increasing  Kic  may  not  be  appropriate.  A number  of  studies  have  reported  that  the 
behavior  of  brittle  and  ductile  materials  is  distinctly  different  and  should  be  addressed 
separately.8^ 

According  to  Lung  and  Mu,  the  variable  Ej  represents  a "step  of  crack  propagation". 
The  authors  do  not  elaborate  on  precisely  how  this  "step"  might  function  in  the  fracture 
process.  In  their  original  work  Lung  and  Mu  suggest  that  fracture  surfaces  are  fractal,  but 
provide  no  model  for  propagation  of  a fractal  pattern  on  a fracture  surface.  No  suitable 
physical  models  were  presented  by  these  or  any  other  authors  at  a later  date. 

The  mathematical  relationship  developed  by  Lung  and  Mu  may  used  to  approximate  a 
value  for  8j.  The  data  shows  that  this  step  length  dimension  is  in  the  order  of  2.5  x 10  7 m 
to  8.0  x 10"5  m. 

Comparison  Against  Mecholsky  et  al's  Model. 

The  model  provided  by  Mecholsky  et  al.  is  the  simplest  to  use.  The  data  obtained 
in  this  study  may  be  readily  modified  to  fit  their  equation  (again,  data  from  set  #1  was  not 
used).  Figure  4-10  shows  the  results  of  this  evaluation.  R2  is  0.77. 
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As  Mecholsky,  Passoja  and  Mackin  argue,  their  formulation  may  be  accepted  on 
the  basis  of  dimensional  analysis.  One  of  the  present  weaknesses  of  fractal  theory  is  that, 
although  fractal  behavior  may  be  identified,  the  relationships  that  are  determined  from 
measurements  (Richardson  graphs  or  dimensional  analysis  graphs)  do  not  identify  the 
actual  underlying  mechanism  which  created  it.  In  developing  their  equation,  Mecholsky 
et  al.  relied  partially  on  dimensional  analysis.  Recent  work  shows  that  their  results  may 
be  well  founded  in  atomic  theory.123 

Equation  4-3  contains  the  term  a<)  which  according  to  the  authors  represents  a 
structure  parameter.  As  noted  previously  in  Figure  2-6,  this  characteristic  length  is 
consistent  with  a free  volume  or  the  sizes  of  glass-crystal  clusters,  stretched  bonds  or 
glass  phase  stretched  bonds.34  The  approximate  dimensions  range  from  1 to  80 
Angstroms.  The  actual  meaning  of  ao  or  its  place  in  the  fracture  process  remains 


unknown. 
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A similar  set  of  assumptions  as  those  made  by  Mecholsky  and  co-workers  were 
employed  here.  That  is,  based  on  the  favorable  results  of  the  Kic  v.  (D*)1 2 relationship, 
it  was  assumed  that  equation  4-3  may  also  represent  the  data  obtained  in  this  report. 
Beginning  with  equation  4-3  and  solving  for  ao  shows  that 

(4-4)  ao  = (K|C2)/(E2  D*) 

A graph  of  Kic  v.  ao' 2 is  illustrated  in  Figure  4-11. 


(ao)12  v K|C 


(ao)12(m1/2) 


Figure  4-11.  Graph  of  ao' 2 v.  Kic. 

(Note:  In  order  to  avoid  confusion,  the  variable  ao  will  be  used  throughout  this 
study  to  represent  the  "structure  parameter",  however,  be  aware  that  ao  may  be  a constant 
or  a variable  depending  on  the  context  in  which  it  is  discussed.). 

According  to  4-3,  ao  increases  with  increasing  K|C.  The  range  of  ao  values  lies 
between  8.2  x 10  s m (at  Kic  = 25.5  MPa  Vm)  and  2.9  x 10‘h  m (at  Kic  = 121.2  MPa  Vm). 
Recall  that  Mecholsky  and  co-workers  have  interpreted  ao  to  be  a structure  dependent 
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dimension  associated  with  the  fracture  process  in  ceramics.  In  their  case  the  values  of  ao 
ranged  between  1 and  80  angstroms.  Considering  that  metallic  alloys  plastically  deform 
at  the  crack  tip  during  loading  and  that  this  process  involves  hundreds  or  perhaps 
thousands  of  atoms,  it  should  not  be  surprising  that  the  values  of  ao  are  in  the  range  of 
tens  of  nanometers  to  microns  rather  than  angstroms.  In  this  case,  ao  may  also  be 
interpreted  as  a structure  parameter.  The  question  remains,  however,  what  feature  is  ao 
representing?  The  continuation  of  this  document  attempts  to  provide  some  insight. 
Modified  Mecholsky  et  al. 

Figure  4- 1 0 showed  a linear  relationship  between  Kic  and  D*.  It  was  suggested 
that  the  relationship  may  be  adequately  described  by  Equation  4-3  or  4-4.  Figure  4-10 
may  also  be  adequately  described  by  a slightly  modified  version  of  Equation  4-3.  In  view 
of  the  fact  that  the  relationship  appears  to  be  linear  and  that  dimensional  requirements 
must  be  satisfied,  Kic  and  D*  may  also  be  related  as  follows 

(4-5)  K|C  = 250  - 402X 

where  the  variables  and  constants  are  interpreted  as  follows: 

Kic  = Plane  strain  fracture  toughness. 

250=  A constant  (MPa  m12). 

402  = A second  constant  composed  of  Eao'  2 (MPa  m1 2). 

E = Young's  modulus. 

ao  = A material  dependent  structure  parameter. 

X = (D*)1  2 = Surface  fractal  dimension. 

In  accordance  with  this  relationship  Eao12  = 402  MPaVm.  Using  the  fact  that  E is 
approximately  200  Gpa  for  heat  treated  steels,  ao  may  be  solved  to  be  equal  to  4.0  x 10"6 


meters. 
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Interpretation  of  a0  and  D* 

The  results  of  this  study  indicate  that  fracture  surfaces  of  metallic  materials  that  fail  in 
a ductile  manner  generate  in  a fractal  pattern.  The  relationships  identified  between 
fracture  toughness  and  D*  suggest  that  the  pattern  is  related  in  some  way  to  the 
characteristics  of  the  plastically  deformed  region  ahead  of  the  crack  tip.  The  processes 
that  control  K)C  take  place  in  this  in  this  zone. 

Three  models  have  been  presented  to  account  for  the  observed  behavior.  In  two  of 
the  models  a structure  parameter,  ao,  may  be  taken  to  be  a variable  while  in  the  third  it  is 
a constant.  Table  4-3  summarizes  some  of  the  principal  findings. 

Attempting  to  identify  which,  if  any,  of  the  models  suggested  is  correct  requires 
knowledge  of  the  exact  mechanism  of  ductile  fracture  which  continues  to  be  investigated. 
The  findings  of  other  researchers  in  the  field  of  fracture  toughness  - microstructure 
relationships,  however,  may  be  useful  in  deciding  what  factors  are  relevant  and  how  they 
are  thought  to  be  related  to  Kic. 

Table  4-3 


Summary  of  Principal  Findings  Derived  from  Mathematical  Models 


Model 

a0(or  s)  range 

K|C-max 

(m) 

(MPaVm) 

Mu  and  Lung 

2.5x1  O' ; to  8.0x1  O'5 

247-257 

Mecholsky  et  al. 

8.2x1  O'8  to  2.9x1  O'6 

250 

Mod.  Mecholsky  et  al. 

4.0x1  O'6 

250 

Characteristic  Distance  Ahead  of  a Crack  Tip. 

It  is  well  known  that  a highly  stressed  field  develops  ahead  of  a crack  tip  upon 
loading.  Calculations  and  computer  simulations  have  shown  that  the  magnitude  and 

distribution  of  the  stresses  may  be  given  by  graphs  such  as  that  shown  in  Figure  4-12,'  24‘ 

1 26 

where  oyy  is  the  stress  approximately  perpendicular  to  the  crack  plane,  o0  is  the  yield 
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stress,  r is  a distance  ahead  of  the  crack  tip,  n is  the  Ramberg-Osgood  strain  hardening 
coefficient  and  J is  the  J integral  measure  of  fracture  toughness. 

The  drop  in  stresses  close  to  the  crack  tip  is  due  to  the  presence  of  a free  surface  and 
significant  crack  tip  blunting.  The  strains  in  the  region  close  to  the  crack  tip  are  also 
considerable.  A typical  plot  for  strains  at  0=  0°  is  illustrated  in  Figure  4-13. 127 
5 
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Some  theories  suggest  that  crack  initiation  occurs  when  a certain  critical  stress  (of)  or 
critical  strain  (ep.crjt)  is  exceeded  over  a characteristic  distance  (Lq*)127'130  ahead  of  the 
crack  tip.  These  models  unite  macroscopically  observed  fracture  toughness  with 
microscopically  observed  material  features  (second  phase  particles,  inclusions,  grain  size 
distribution  etc.).  In  brittle  materials  the  fracture  criterion  is  based  on  the  critical  stress 
while  in  ductile  materials  the  criterion  is  founded  on  reaching  a critical  strain.  When 
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Figure  4-12.  Stress  Field  Ahead  of  a Crack  Tip  [Ref  124], 

Copyright  (1995)  from  Fracture  Mechanics  by  T.L.  Anderson.  Reproduced  by 
permission  of  Routledge/Taylor  & Francis  Group,  LLC. 
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each  of  these  critical  quantities  is  satisfied,  over  a characteristic  distance,  the  crack  front 
advances. 

A number  of  concepts  and  models  involving  a characteristic  distance  have  been 
suggested  for  brittle  fracture.  Considering  that  the  subject  of  this  work  involves  primarily 
ductile  fracture  the  reader  is  referred  to  other  studies  for  the  details  of  these 
approaches.1'7’  I‘9, l30'133  For  the  purposes  of  this  work  it  is  sufficient  to  recognize  that 
fracture  is  predicted  to  occur  in  brittle  alloys  when  a given  critical  stress  is  reached  over  a 
characteristic  distance. 


Figure  4-13.  Strain  Field  Ahead  of  a Stressed  Crack  Tip  [Ref.  127], 

More  useful  within  the  context  of  this  study  is  the  fracture  toughness-microstructure 
relationship  suggested  for  ductile  fracture  conditions.  In  this  case,  the  dominant 
fractographic  features  are  microvoids.  A number  of  models  have  surfaced  which  appear 
to  have  predictive  abilities.  Of  these,  one  of  the  more  successful  ones  appears  to  be  the 
stress-modified  critical  strain  (SMCS)  criterion.127, 128, 131, 134, 135 
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According  to  the  SMSC  fracture  takes  place  when  the  critical  strain  (sp.cnt)  reaches  a 
critical  value  over  a critical  distance  (L0*).  The  value  of  sp.crit  is  given  by  equation  4-6. 
The  value  or  meaning  of  Lo  is  defined  differently  by  various  researchers. 

(4-6)  8pcntlcal  = a exp(P  am/ae) 

where 

^critical  _ crjtjcai  strajn  ahead  of  the  crack  tip  at  failure. 

a and  p = Experimentally  determined  parameters. 

am  = Hydrostatic  stress  (a.k.a.  mean  stress). 

ae  = Von  Mises  stress  (a.k.a  effective  stress). 

Note:  For  many  steels  P is  taken  to  be  -1.5. 

These  concepts  were  examined  extensively  by  Ritchie  and  co-workers.128  They 
incorporated  the  SMSC  into  a equation  relating  fracture  toughness  to  the  critical  strain 
and  critical  distance  as  follows 


(4-7)  K,c  = C(E'oo£fV)1/2 


where 

Kic  = Plane  strain  fracture  toughness. 

C = A constant. 

Sf*  = Critical  strain  ahead  of  the  crack  tip  at  failure. 

Go  = Flow  stress. 

E'  = Young's  Modulus  under  small  scale  yielding  conditions, 
and 


(4-8)  sf*  = ln(dp/Dp) 

0.28  exp(  1 .5  am/Ge) 


where 

£f*  = Critical  fracture  strain. 

om  = Hydrostatic  stress  (a.k.a.  mean  stress). 

oe  = Von  Mises  stress  (a.k.a  effective  stress). 

dp  = Mean  spacing  between  void  nucleating  inclusions. 

Dp  = Void  initiating  particle  diameter. 


One  of  the  difficulties  with  this  type  of  model  (equation  4-7)  is  that  Lo*  is  often  not 
known  or  not  readily  obtainable  from  either  laboratory  tests  or  metallographic 
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measurements.  Lo*  is  usually  selected  on  the  basis  of  curve  fitting  or  computer  based 


numerical  modeling. 


In  a study  on  high  strength  nuclear  pressure  vessel  steels  (ASME  SA533B-1  in  a 


bainitic  heat  treat  condition.),  Ritchie  et  al.l2fi  estimated  L0*  to  be  in  the  order  of  50  to  100 


pm  for  low  toughness  steels  (KiC  in  the  order  of  40  Mpa  V m)  and  300  to  350  pm  for  high 


toughness  steels  (Kic  in  the  order  of  225  MPaVm.).  Their  findings  are  shown  in  Figure  4- 


14.  A brittle  alloy  toughness  model  (also  known  as  RKR  theory)  was  used  for  the  lower 


shelf  curves  while  SMCS  was  applied  to  the  upper  shelf  data. 
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Figure  4-14.  Fracture  Toughness  v Temperature  for  ASME  SA533B-1, 

Cl  2,  bainitic  condition.  Brittle  and  ductile  models  of  fracture  initiation 
showing  L0*  fitted  to  the  data.  {Ref  128], 

While  the  practical  difficulties  of  obtaining  Lo*  and  Sf*  limit  the  use  of  equations  4-6 


to  4-8  in  real  engineering  applications,  it  does  help  to  reason  through  trends  in  fracture 


toughness  due  to  variations  in  stress  state  (triaxiality  = am/Ge),  material  characteristics 


(flow  stress  - oo)  and  material  "cleanliness"  (reflected  in  ratio  of  dp/Dp).  The  SMCS 
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model  has  been  used  in  conjunction  with  FEA  analyses  to  numerically  predict  changes  in 
fracture  toughness  as  a result  of  welding,  forming  or  plastic  deformation  induced  residual 
stresses.136' 137 

In  a recent  report  Kanvinde  and  Deierlein  have  incorporated  SMCS  theory  into  a 
numerical  model  for  predicting  ductile  fracture  in  several  structural  steels.138’139  As  in 
previous  uses  of  this  model  they  calculate  the  critical  strain  using  equation  4-6.  In  a 
unique  manner,  however,  they  infer  Lo  from  fractographic  images  rather  than  from  curve 
fitting.  Their  model  assumes  that  the  critical  step  for  unstable  fracture  occurs  when 
several  large  holes  (void  clusters)  formed  by  the  nucleation  and  growth  of  smaller 
dimples  link-up  thereby  initiating  fast  fracture.  The  diameter  of  the  large  void  clusters  in 
the  fracture  plane  is  taken  to  represent  Lo*. 

Despite  the  different  methodology  and  steel  alloys,  the  results  of  Kanvinde  and 
Deierlein1 17  are  comparable  with  those  of  Ritchie  et  al.128  The  critical  lengths  in  the 
upper  shelf  are  100  to  300  micro-meters  and  300  to  350  micrometers,  respectively. 

As  illustrated  by  the  work  of  the  two  previous  research  groups,  the  selection  and 
meaning  of  Lo  is  not  clear.  Other  researchers  have  suggested  additional  interpretations 
of  the  critical  distance  in  ductile  materials.  Panontin  et  al.134  and  Ritchie  et  al. 128,130  have 
suggested  that  the  length  scale  is  related  to  grain  size  (usually  2 to  4 times  the  mean  grain 
diameter).  Alternatively,  Chi  et  al.140  and  Rousselier141  suggest  that  the  critical 
dimension  is  attributable  to  microvoid  sizes  and  distance  between  inclusions. 
Relationship  Between  L0  , D*  and  a0. 

Recall  that,  depending  on  the  model,  the  results  of  the  present  study  calculated  a<> 
values  in  the  order  of  0.08  pm  to  80  pm  (Table  4-3).  These  values  are  anywhere  from  2 
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to  3,500  times  smaller  than  those  proposed  by  Ritchie  et  al.128  and  Kanvinde  and 
Deierlein. 1 38 

Only  one  of  the  models  (Mu  and  Lung;  s ~ 8.0x1 0°)  suggests  a step  size  for  crack 
growth  which  is  comparable  to  the  reported  range  of  characteristic  distances  (Lo* 
between  100  to  350  microns).  The  remaining  models  indicate  that  ao  is  considerably 
smaller  that  Lo*.  As  a result,  these  two  quantities  are  considered  to  represent  different 
but  related  physical  distances. 

Considering  the  structure  parameter  in  a fractal  context  leads  to  the  concept,  as  other 
researchers  have  reasoned,74'83'87  that  a « may  represent  or  is  related  to  the  step  or  steps 
involved  in  the  generation  of  a fractal  surface.  The  variable  Lo*,  on  the  other  hand, 
represents  the  final  size  of  the  microstructural  damage  achieved  just  prior  to  fracture. 
Ritchie  et  al.  have  noted  that,  in  ductile  materials,  the  joining  of  several  dimples  in  the 
crack  tip  plastic  zone  is  required  prior  to  crack  front  advancement.128 

The  data  in  this  study  has  been  found  to  be  represented  by  equation  4-3  as  follows 

KIC  = E(aoD*),/2 

Note  that  this  function  is  similar  to  that  discussed  by  Pandey  and  Banerjee,138  and 
Ritchie  et  al.1'8  for  ductile  fracture 

(4-9)  KIC  = ((Constant)E'ayid£f*L0*)1 2 

Observing  that  the  yield  strength,  oyid,  is  equal  to  Eeyid,  E'  = E/(l-v2)  for  plane  strain 
and  incorporating  1/(1-  v2)  into  the  constant,  then 

(4-10)  E(aoD*)1/2  = E(C£yld£fV)1/2 

and 


(4-11)  a<)D*  = C8yid£f*L0* 
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or 


(4- 12a)  ao  = (CsyidEf  Lo*)/D*  and  (4- 12b)  D*  = (Ceyid8f*Lo*)/ao 
Equations  4- 12a  and  4- 12b  suggest  that  there  is  a close  interrelationship  between 
parameters  representing  material  properties  (eyid,  £f*)  and  those  representing  structures 
involved  in  the  dimple  rupture  damage  evolution  process  (ao,  D*  and  Lo*). 

Relationship  Between  D*,  ao  and  Microstructural  features. 

Earlier  in  this  report  it  was  observed  that  Garrison  et  al.  have  performed  a number 
of  studies  showing  the  relationship  between  fracture  toughness  and  the  product  of  the 
particle  spacing  - void  growth  ratio  (X0(RV/R,))"  M 16  (see  Fig.  4-8).  Garrison  has  also 
considered  the  possible  relationship  between  his  findings  and  those  expressed  in  equation 
4-9.  He  has  proposed  the  following  integration  of  the  two  concepts' 1 1 
Replacing  Sf*L0  with  X0(RV/Ri) 

(4-13)  Kic  = ((Constant)Xo(Rv/Ri)Eayid)12 
In  addition,  he  proposes  that  one  may  reasonably  go  one  step  further  and  state  that 


if 


(4-14)  8f  Lo  - aXo(Rv/Rj)|R0 


then 


(4-15)  ef*  = b { ( Rv/Rj |Ro) } “ 

and 


(4-16)  Lo*  = cX0{(Rv/Rj)|Ro}(1'o) 

Where  the  constants  are:  a = cb  and  c > 1 

Garrison  notes  that  these  equalities  are  consistent  with  the  physical  meaning 
previously  ascribed  to  these  variables.1 1 1 For  instance,  s/  has  been  considered  to 
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designate  a fracture  strain  which  may  be  correlated  to  a measure  of  the  void  growth  ratio. 
The  characteristic  distance,  Lo  , in  turn,  has  been  suggested  to  designate  a distance 
between  fracture  initiating  particles  or  grain  boundaries  which  may  also  be  related  to  the 
inclusion  spacing  X0. 

Following  the  same  reasoning  suggested  by  Garrison,  the  findings  of  this  report  may 
be  equated  to  microstructural  features  for  a limited  range  of  applicability1 11,115 

That  is 

(4-17)  KIC  = E(aoD*)l/2  = E(Cey,defV)l/2  = E(CaX0(Rv/Ri|Ro)syid)l/2 

Where  C and  Ca  are  constants. 

Considering  the  findings  of  previous  research,  fracture  toughness  appears  to  be 
closely  related  to  two  major  microstructural  features;  the  distance  between  void 
nucleating  particles  and  the  degree  of  matrix  ductility.  The  model  which  includes  D*  and 
best  describes  these  features  is  that  suggested  by  Mecholsky  et  al.30  (equation  4-3).  Note 
that  Mecholsky's  group  considered  ay  to  remain  constant  for  a material  class.  While  that 
may  be  true  for  ceramics  it  does  not  appear  to  be  so  for  metals.  Recognizing  that  heat 
treatments  will  change  particle  sizes  and  distances,  as  well  as  matrix  properties,  a two 
parameter  description  involving  material  and  surface  features  is  thought  to  be  necessary 
to  determine  fracture  toughness.  A single  parameter  description  for  K|C  based  solely  on 
the  surface  geometry  does  not  seem  reasonable  at  this  time. 

Table  4-4  shows  the  approximate  range  of  dimensions  for  ao  against  the  specimen 
numbers.  Note  that  the  more  brittle  specimen  (low  numbers)  exhibit  small  values  of  ao 
and  that  these  values  increase  with  toughness.  These  dimensions  are  consistent  with  the 
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range  of  sizes  expected  in  the  average  voids,  particle  sizes  and/or  interparticle  spacing  of 
second  phases  in  AISI  4340  heat  treated  to  various  strength  levels. 

Table  4-4 


Specimen  Number  v.  a0  Dimensions 


Specimen  Numbers 

ao  (meters) 

1-6 

-5  to  9.5  x 

7-15 

-2  to  9.5  x 10~7 

16-24 

~1.3  to  3.0  x 10'6 

Explaining  Fractal  Behavior 

As  observed  in  equation  4-17  the  mathematical  relationships  presently  available  to 
define  D*  and  ao  do  not  adequately  explain  the  precise  meaning  of  these  variables, 
however,  they  do  suggest  that  they  are  related  to  each  other  and  do  appear  to  represent 
microstructural  characteristics.  The  structure  parameter  appears  to  be  closely  related  to 
some  significant  structure  dimension  while  the  behavior  of  D*  is  consistent  with  the 
alloy's  flow  behavior. 

It  has  been  discussed  that  AISI  4340  as  well  as  many  other  high  strength  alloy  steels 
fracture  by  a process  of  microvoid  nucleation,  growth  and  coalescence.  Nucleation 
begins  at  the  largest  particles  (MnS  and  other  inclusions).  As  additional  strain  is  applied 
smaller  particles  (cementite  and  other  carbides)  begin  to  participate  in  the  nucleation 
process.  Growth  of  the  initial  voids  leads  to  localization  of  strains  between  cavities.  The 
strained  regions  nucleate  additional  voids  (a  void  sheet)  which  eventually  connects  the 
larger  voids. 

A completely  separated  fracture  surface  reveals  the  presence  of  dimples  of  ever 
decreasing  size.  Ishikawa14-  has  suggested  that  this  idealized  pattern  of  circles  adjacent 
to  smaller  circles  is  fractal  in  nature  (see  Fig  4-15).  The  results  of  this  study  agree. 
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Ishikawa,  however,  as  most  researchers  in  this  area,  was  not  able  to  explain  the  roots  of 
this  behavior . 

Fracture  surface  observations  made  throughout  this  research,  the  data  collected,  and 
the  findings  of  previous  researchers  suggests  that  fracture  surfaces  are  created  by  the 
application  of  some  prescribed  set  of  rules.  One  of  the  rules  seems  to  be  that  dimples 
should  have  nearly  the  same  proportions  within  a given  range  of  void  nucleating  particle 
sizes. 


Figure  4-15.  Ishikawa's  idealized  view  of  a 
ductile  fracture  surface  dimple  pattern. 

Observe  dimples  adjacent  to  dimples  of  ever 
decreasing  size.  [Ref.  142] 

A fractal  surface  requires  that  the  observed  fractal  pattern  appear  similar  at  a wide 
range  of  magnifications.  In  other  words,  two  images  of  dimple  rupture  taken  at  different 
magnifications  from  the  same  surface  should  appear  the  same  or  nearly  so.  This 
requirement,  in  turn,  implies  that  void  sizes  should  follow  some  degree  of  proportionality 
from  largest  to  smallest. 

As  noted  previously.  Garrison  et  al.  have  examined  a number  of  high  strength  steel 
fracture  surfaces.  Their  findings  show  that  one  of  the  features  that  appears  to  remain 
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relatively  constant  throughout  a range  of  inclusion  sizes,  as  well  as  being  an  indicator  of 
fracture  toughness,  is  the  ratio  of  the  void  radius  to  the  inclusion  radius  (Rv/R,).  Several 
of  their  graphs  are  provided  in  Figures  4- 1 6 to  4- 1 9.  The  data  shows  that  for  the  same 
range  of  features  (dimple  sizes)  sampled  by  the  slit  islands  used  in  this  study  (~0.8 
microns  and  larger),  the  void  diameters  were  nearly  constant.  Although  the  steels 
presented  in  these  figures  are  different  than  AISI  4340,  the  microvoid  coalescence 
process  is  though  to  be  controlled  by  the  same  basic  factors  (inclusion  sizes/distribution 
and  matrix  flow  properties)  and  develops  in  the  same  manner  (nucleation,  growth  and 
coalescence  by  void  sheets). 


Figure  4- 16.  Graph 
illustrating  the  void 
growth  ratio  (Rv/R,) 
as  a function  of  the 
inclusion  radius  for  a 
high  strength  Ni-Si- 
MoV  steel.  [Ref. 
1121. 


Figure  4- 17.  Graph 
illustrating  the  void 
growth  ratio  (Rv/R,) 
as  a function  of  the 
inclusion  radius  for  a 
high  strength  Ni-Si 
steel.  [Ref.  112] 
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Figure  4- 18.  Graph 
illustrating  the  void 
growth  ratio  (Rv/Ri) 
as  a function  of  the 
inclusion  radius  for  a 
high  strength  AF 
1410  steel  aged  at 
425°C.  [Ref.  113] 
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Figure  4-19.  Graph 
illustrating  the  void 
growth  ratio  (Rv/R,) 
as  a function  of  the 
inclusion  radius  for  a 
high  strength  AF 
1410  steel  aged  at 
5 10°C1.  [Ref.  113], 
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Additional  evidence  that  particle  sizes  and  the  associated  void  growth  proportions 
may  account  for  fractal  behavior  may  be  found  in  Figure  2-13.  This  graph  shows  that  for 
the  range  of  particle  sizes  thought  to  be  most  responsible  for  toughness  (1  to  100  micron 
range),  the  fractal  dimension  is  constant.  Again,  it  is  noted  that  the  steels  evaluated  by 
Garrison1 1-1 13  and  Dauskardt  et  al.83  are  different  but  the  underlying  fracture  mechanism 
is  thought  to  be  the  same. 

Figure  2-13  also  shows  that  a different  fractal  dimension  (lower  value)  is  observed  at 
very  small  step-lengths.  At  these  smaller  scales  the  factors  that  are  influencing  the 
fracture  surface  topography  have  been  attributed  to  slip.  Note  that  at  the  small  particle 
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sizes  Garrison  et  al.  found  an  increase  in  the  void  growth  ratio.  That  is,  as  suggested  by 
this  study,  a more  stretched  out  surface  results  in  a decreased  fractal  dimension. 

Changes  in  D*  (surface  topography)  With  Microstructure. 

The  fractal  dimension  of  a fracture  surface,  as  well  as  fracture  resistance,  appears  to 
be  closely  related  to  the  flow  properties  of  a steel  alloy.  Microstructural  changes  which 
affect  toughness  also  change  the  fracture  surface  topography.  As  noted  in  the  previous 
sections,  D*  seems  to  be  tied  to  strain  to  yield  as  well  as  the  plastic  fracture  strain. 

Plastic  fracture  strain  is  reflected  in  the  degree  of  void  growth  from  second  phases.  Other 
researchers  agree  that  the  extent  of  void  deformation  provides  a reliable  indication  of 
alloy  toughness.143, 144 

It  has  been  reported  that  higher  toughness  alloys  display  deeper  and  wider 
microvoids.14"  This  study  has  again  confirmed  that  the  fracture  surfaces  of  tougher  alloys 
shift  from  a series  of  jagged  interconnected  planes  (brittle  behavior)  to  a series  of 
relatively  smooth,  parabolic,  cavities  (ductile  behavior).  It  is  suggested  that  this  shift  is 
responsible  for  the  observed  changes  in  D*.  Lower  values  of  D*  are  taken  to  indicate  a 
smoother  more  stretched  out  surface.  Although  on  a macroscopic  level  the  surface  is 
rougher,  displaying  greater  peaks  and  valleys,  on  a microscopic  level  the  surface  is 
becoming  smoother.  This  is  consistent  with  a large  number  of  studies  characterizing  the 
relationship  between  toughness  and  dimple  geometry.  Particularly,  those  involving 
steel. 104-1 16145  Figures  4-20  to  4-33  illustrate  the  more  elongated  surface  appearance  of  a 
high  toughness  heat  treatment  condition  compared  against  that  of  two  lower  toughness. 


less  ductile,  conditions. 
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Figure  4-20.  Fracture  surface  from  specimen  #20  (42  HRC)  at  -200X. 
Many  dimples  noted. 
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Figure  4-24.  Fracture  surface  #20  (42  HRC)  at  ~20,000X. 


Figure  4-25.  Fracture  surface  #20  (42  HRC)  at  ~50,000X. 
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Figure  4-26.  Fracture  surface  from  specimen  #7  (49.7  HRC)  at  -200X. 
Many  dimples/cleavage  noted. 


Figure  4-27.  Fracture  surface  #7  (49.7  HRC)  at  ~1,000X. 
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Figure  4-32.  Image  of  specimen  # 6.  Shallow  dimples  are  beginning  to 
form  on  the  fracture  surface. 


Figure  4-33.  Image  of  specimen  # 6.  Shallow  dimples  are  beginning  to 
form  on  the  fracture  surface. 


CHAPTER  5 
CONCLUSIONS 

The  results  of  this  study  provide  the  basis  for  forming  the  following  conclusions: 

1 ) The  fracture  surfaces  of  AISI  4340  steel  heat  treated  to  achieve  various  values  of 
fracture  toughness  exhibit  fractal  qualities.  When  measured  using  the  Richardson  slit- 
island  method,  the  fracture  surfaces  approximate  the  mathematical  characteristics  of  a 
fractal  surface.  In  particular,  they  provide  linear  relationships  in  the  log( yardstick)  v. 
log( total  length)  graphs  and  a fractional  dimensional  increment  in  the  expected  range  of 
values  (0  < D*  < 1). 

2)  Fracture  surfaces  exhibit  statistically  fractal  qualities.  Fracture  surfaces  are  formed 
from  real  materials  containing  a distribution  of  defects,  grain  sizes,  grain  orientations, 
micro-chemistries,  inclusions  and  second  phases.  Therefore,  they  do  not  form  exact 
fractal  structures  (not  comparable  to  mathematical  fractals). 

3)  Fracture  surfaces  are  not  isotropic  or  homogeneous,  therefore,  many  slit-island 
measurements  are  necessary  to  obtain  statistically  significant  values  of  D*. 

4)  The  fracture  surfaces  in  this  study  exhibited  a decrease  in  fractal  dimension  with  an 
increase  in  fracture  toughness  for  materials  that  failed  in  a ductile  manner.  The  results 
are  consistent  with  the  results  of  the  majority  of  similar  studies. 

The  fracture  surfaces  that  showed  primarily  cleavage  topography,  corresponding  to 
more  brittle  type  fracture,  were  consistent  with  an  increase  in  the  fractal  dimensional 
increment  with  an  increase  in  fracture  toughness,  as  observed  for  other  brittle  materials. 
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5)  The  decreasing  values  of  D*  with  increasing  fracture  toughness  may  be  modeled  by 
equations  provided  by  Mu  and  Lung  or  Mecholsky,  Passoja  and  Mackin.  Both  of  these 
functions  suggest  the  involvement  of  a structure  parameter  (ao  or  s),  having  units  of 
length,  on  the  fracture  process.  The  equation  provided  by  the  latter  group,  however, 
seems  to  be  more  consistent  with  other  research  on  ductile  fracture. 

6)  The  significance  of  the  structure  parameter,  ao,  for  materials  that  fail  in  a ductile 
manner  is  suggested  by  a fractal  surface  generation  model.  This  model  proposes  that 
fractal  fracture  surfaces  are  produced  by  the  nucleation,  growth,  and  coalescing  of 
dimples  of  ever  decreasing  size.  The  sequence  is  that  of  large  dimples  nucleating  at  the 
largest  inclusions  followed  by  the  continuous  growth  of  ever  decreasing  dimples  adjacent 
to  them.  In  this  context,  ao,  represents  a mean  interparticle  spacing,  a mean  void  diameter 
or  similar  structural  parameter. 

7)  In  the  proposed  fractal  fracture  surface  generation  model,  D*  is  taken  to  represent  the 
degree  of  ductility  available  in  the  alloy  matrix  (the  material  that  surrounds  the  primary 
voids).  An  increase  in  fracture  toughness  is  found  when  dimples  are  more  stretched  out. 
This  surface  geometry  seems  to  be  consistent  with  a lower  value  of  D*. 

8)  SEM  fractographs  of  (ductile)  low  fracture  toughness  and  high  fracture  toughness 
specimens  exhibit  dimples  and  patterns  consistent  with  the  proposed  descriptions  for  ao 
and  D*. 

9)  The  model  is  consistent  with  other  models  for  ductile  fracture  which  state  that  a 
critical  strain  must  be  reached  over  a critical  distance  in  order  for  fracture  to  proceed. 

The  fractal  model  involves  the  initial  void  nucleation,  growth  and  coalescence  phases  just 
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prior  to  instability.  When  the  pattern  reaches  a critical  size,  the  crack  front  advances  in  a 
stable  or  unstable  manner. 

10)  The  model  is  able  to  qualitatively  explain  trends  in  fracture  toughness  and  expected 
behavior  when  microstructure  and  composition  are  modified.  Factors  which  increase  the 
product  of  ao  and  D*  are  expected  to  increase  KiC.  K.1C  rises  when  the  initial  fractal 
generator  step  sizes  (dimple  sizes  or  dimple  array  size)  increase  and  the  surface  is 
stretched  out.  Both  of  these  requirements  results  in  very  high  energy  consumption  in  the 
formation  of  a fracture  surface. 


CHAPTER  6 
FUTURE  WORK 


1 ) This  study  provides  some  possible  answers  to  the  question  of  whether  a fracture 
surface  is  fractal,  why  such  a surface  may  become  that  way  and  what  relationship  this 
may  have  to  fracture  toughness.  It  has  been  suggested  that  KiC,  ao  and  D*  may  be  related, 
at  least  in  part,  to  interparticle  spacing  and  void  growth  proportions.  These  suggestions 
should  be  followed  up  with  feature  measurements  which  either  confirm,  modify  or 
contradict  the  results  of  this  study. 

2)  This  work  suggests  that  Kic  = E(aoD*)' 2 does  represent  the  data  accurately.  More 
work  with  other  alloy  systems  is  required  to  determine  if  this  relationship  might  be 
universal  or  only  appropriate  for  a single  alloy  class. 

3)  The  results  of  this  study  may  be  useful  in  determining  how  best  to  construct  a model 
for  how  fracture  surfaces  develop.  It  might  be  useful  to  do  so  with  interparticle  spacings 
and  rules  governing  void  growth  proportions  as  part  of  the  input  data. 

4)  The  fractal  dimension  was  determined  for  features  in  the  range  of  1 to  100  microns.  It 
would  be  of  interest  to  determine  how  D*  changes  at  lower  scales.  This  work  and  that  of 
other  researchers  suggests  that  D*  will  decrease  if  the  features  at  smaller  scales  are 
smoother  or  more  stretched  out  than  those  at  the  greater  scales. 
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Load  Ibf 


APPENDIX  A 

Examples  Of  Load  Displacement  Diagrams  Used  In  Fracture  Toughness  Measurements 


Figure  A- 1.  Load 
displacement 
diagram  for 
specimen  #2. 

Kic  = 27.6  MPaVm. 


Figure  A-2.  Load 
displacement 
diagram  for 
specimen  #12. 

Kic  = 57.6  MPaVm. 
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Figure  A-3.  Load 
displacement  diagram 
for  specimen  #16. 

Kq  = 99.3  MPaVm 


Figure  A-4.  Load 
displacement  diagram 
for  specimen  #17. 

Kq  = 96.4  MPaVm 


Load  Ibf  Load  Ibf 
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Figure  A-5.  Load 
displacement 
diagram  for 
specimen  #18. 

Kq  = 85.9  MPaVm 


Figure  A-6.  Load 
displacement 
diagram  for 
specimen  #19. 

Kq=  116.3  MPaVm 


Load  Ibf  Load  Ibf 
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Figure  A-7.  Load 
displacement 
diagram  for 
specimen  #20. 

K0=  1 18.4  MPaVm. 


Figure  A-8.  Load 
displacement 
diagram  for 
specimen  #21. 

Kq=  105.9  MPaVm. 


Load  Ibf 
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Figure  A-9.  Load 
displacement 
diagram  for 
specimen  #22. 

K0=  120.5  MPaVm 


Figure  A- 10.  Load 
displacement  diagram 
for  specimen  #23. 
Kq=  121.3  MPaVm 


Load  Ibf 
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Figure  A- 1 1 . Load 
displacement  diagram 
for  specimen  #24. 
Kq=  1 15.5  MPaVm 


APPENDIX  B 

Textbook  Images  Used  To  Verify  Fractal  Measurements 


Figure  B-l . Computer  generated  fractal  surfaces  used  to  verify 
performance  of  fractal  image  measurement  software. 


(a) 


(b) 


Figure  B-2.  (a)  Classic  Koch  triadic  (D*  - 
[Images  modified  slightly  to  fit  this  page] 


.262)  and  (b)  quadratic  (D*  = .500)  curves. 


Reference  for  all  three: 

J.C.  Russ  "Fractal  Surfaces"  (Plenum  Press,  New  York,  1994)  pg.  4,  5 and  31. 
(Reproduced  with  kind  permission  of  Springer  Science  & Business  Media) 
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APPENDIX  C 

Slit  Island  Images  Used  In  Fractal  Analysis 


The  images  contained  within  this  Appendix  illustrate  the  slit-island  geometries 
obtained  from  the  fracture  toughness  specimens  and  employed  in  fractal  dimension 


measurement.  The  images  are  organized  in  order  of  lowest  to  highest  fracture  toughness. 


The  data  reported  in  Table  4-1  is  reproduced  again  here  for  ease  of  reference. 


Table  C-l.  Fracture  Toughness  - Fractal  Dimension  Results  (Raw  Data) 


Sample 

Number 

Hardness 

Before 

(HRC) 

Hardness 

After 

(HRC) 

Crack 

Length 

m 

Fracture 

Load 

Newtons 

Shape 

Factor 

Fracture 
Toughness 
Mpa  V m 

D* 

1 

56.3 

56.3 

0.02667 

6978.9 

10.45 

25.6 

0.199 

2 

56.4 

56.4 

0.02540 

8153.2 

9.66 

27.6 

0.204 

3 

55.7 

55.7 

0.02692 

71 16.8 

10.63 

26.5 

0.187 

4 

55.4 

52.4 

0.02540 

7561.6 

9.66 

25.6 

0.282 

5 

55.5 

52.0 

0.02464 

8495.7 

9.23 

27.5 

0.264 

6 

55.7 

51.7 

0.02464 

9340.8 

9.23 

30.2 

0.294 

7 

53.9 

49.7 

0.02438 

16680.0 

9.09 

53.2 

0.250 

8 

55.8 

49.6 

0.02540 

14233.6 

9.66 

48.2 

0.237 

9 

56.1 

49.9 

0.02286 

21350.4 

8.34 

62.4 

0.257 

10 

52.5 

46.8 

0.02540 

16457.6 

9.66 

55.7 

0.193 

11 

51.5 

48.0 

0.02616 

16457.6 

10.12 

58.4 

0.219 

12 

52.8 

47.5 

0.02464 

17792.0 

9.23 

57.6 

0.210 

13 

52.5 

45.7 

0.02286 

28467.2 

8.34 

83.3 

0.177 

14 

52.3 

46.5 

0.02426 

21795.2 

8.98 

68.6 

0.221 

15 

51.5 

46.8 

0.02426 

21128.0 

8.98 

66.5 

0.187 

16 

51.8 

43.0 

0.02464 

30691.2 

9.23 

99.3 

0.157 

17 

52.5 

43.8 

0.02540 

28467.2 

9.66 

96.4 

0.121 

18 

53.1 

44.2 

0.02718 

22684.8 

10.80 

85.9 

0.133 

19 

52.6 

41.3 

0.0245 1 

36473.6 

9.09 

116.3 

0.116 

20 

52.0 

42.0 

0.02489 

36028.8 

9.37 

118.4 

0.120 

21 

51.4 

42.1 

0.0273 1 

27724.4 

10.89 

105.9 

0.091 

22 

55.5 

40.2 

0.02540 

35584.0 

9.66 

120.5 

0.176 

23 

50.6 

40.5 

0.02540 

35806.4 

9.66 

121.3 

0.126 

24 

52.3 

40.4 

0.02464 

35673.0 

9.23 

115.5 

0.159 
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Figure  C-l.  A graph  of  K[cv.  D*  illustrating  the  degree  of  variability  in  D*  present  in 
the  data. 
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Figure  C-la 
D*  = 0.20918 


Figure  C-lb 
D*  = 0.1 8508 


Figure  C-lc 
D*  = 0.20299 
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Figure  C-2a 
D*  =0.21019 


Figure  C-2b 
D*  = 0.21092 


Ill 


Figure  C-3a 
D*  = 0.1 7086 


Figure  C-3b 
D*  =0.16807 


Figure  C-3c 
0*  = 0.221 11 
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Figure  C-4a 
D*  =0.31828 


Figure  C-4b 
D*  = 0.28206 


Figure  C-4c 
D*  = 0.24483 
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Figure  C-5a 
D*  =0.29884 


Figure  C-5b 
D*  =0.27498 
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Figure  C-6a 
D*  =0.26536 


Figure  C-6b 
D*  = 0.29401 


Figure  C-6c 
D*  = 0.32258 
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Figure  C-7a  and 
b 

D*(a)  = 0.24242 
D*(b)  = 0.22037 


Figure  C-7c 
D*  = 0.30174 


Figure  C-7d 
D*  = 0.23678 
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Figure  C-8a 
D*  = 0.27602 


Figure  C-8b 
D*  = 0.22414 


Figure  C-8c 
D*  = 0.22633 
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Figure  C-9a 
D*  = 0.33127 


Figure  C-9b 
D*  = 0.23291 


Figure  C-9c  and 
d 

D*(c)  = 0.21328 
D*(d)  = 0.24006 
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Figure  C-9e 
D*  = 0.28716 


Figure  C-9f 
D*  = 0.23661 
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Figure  C-lOa 
D*  = 0.21310 


Figure  C-lOb 
D*  =0.19569 


Figure  C-lOc 
D*  = 0.17148 
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Figure  C-l  la 
D*  = 0.27696 


Figure  C-l  lb 
D*  = 0.20276 


Figure  C-l lc 
D*  = 0.17871 
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Figure  C-12a 
D*  = 0.21841 


Figure  C-12b 
D*  = 0.21836 


Figure  C-12c 
D*  = 0.19293 
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Figure  C-13a 
D*  = 0.20662 


Figure  C-13b 
D*  =0.16948 


Figure  C-13c 
D*  = 0.15626 
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Figure  C-14a 
D*  = 0.22060 


Figure  C-14b 
D*  = 0.2151 1 


Figure  C-14c 
D*  = 0.22641 
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Figure  C-15a 
D*  = 0.21770 


Figure  C-15b 
D*  = 0.15615 


Figure  C-15c 
D*  = 0.18659 
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Figure  C-16b 
D*  = 0.12210 
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Figure  C-17a 
D*  = 0.125 14 


Figure  C-l  7b 
D*  =0.10469 


Figure  C-l  7c 
D*  = 0.13395 
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Figure  C-18a 
D*  = 0.1 1616 


Figure  C-18b 
D*  = 0.14686 


Figure  C-18c 
D*  =0.13636 
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Figure  C-19a 
D*  = 0.13191 


Figure  C-19b 
D*  = 0.14458 


Figure  C-19c 
D*  = 0.07253 
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Figure  C-20a 
D*  = 0.12714 


Figure  C-20c 
D*  =0.10179 
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Figure  C-21a 
D*  = 0.12966 


Figure  C-21b 
and  c 

D*(b)  = 0.08658 
D*(c)  = 0.07294 


Figure  C-21d 
D*  = 0. 1 1991 
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Figure  C-21e 
D*  = 0.04018 


Figure  C-21f 
D*  = 0.09440 
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Figure  C-22a 
D*  = 0.1 5203 


Figure  C-22b 
D*  = 0.18989 


Figure  C-22c 
D*  = 0.18587 


134 


Figure  C-23a 
D*  = 0.09862 


Figure  C-23b 
D*  = 0.09889 


Figure  C-23c 
D*  = 0.1 5203 
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Figure  C-24a 
D*  =0.13333 


Figure  C-24b 
D*  = 0. 14525 


Figure  C-24c 
D*  = 0.19721 
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